AFWL-TR-76-1 67 


AFWL-TR- 

76-167 


♦* 


00 


00 


UNSTABLE  RESONATOR  THEORY 


University  of  Arizona 
Optical  Sciences  Center 
Tucson,  Arizona  85721 


March  1977 


Final  Report 


Approved  for  public  release;  distribution  unlimited. 


AIR  FORCE  WEAPONS  LABORATORY 
Air  Force  Systems  Command 
Kirtland  Air  Force  Base,  NM  87117 


AFV/L-TR-  76-167 


This  final  report  was  prepared  by  the  Optical  Sciences  Center,  University 
of  Arizona,  Tucson,  Arizona,  under  Contract  F29601-75-C-0088,  Job  Order  317J- 
0199  with  the  Air  Force  Weapons  Laboratory,  Kirtland  Air  Force  Base,  New  Mexico 
Dr.  P.  V.  Avizonis  (AR)  was  the  Laboratory  Project  Officer-in-Charae. 

When  US  Government  drawings,  specifications,  or  other  data  are  used  for  any 
purpose  other  than  a definitely  related  Government  procurement  operation,  the 
Government  thereby  incurs  no  responsibility  nor  any  obligation  whatsoever,  and 
the  fact  that  the  Government  may  have  formulated,  furnished,  or  in  any  way 
supplied  the  said  drawings,  specifications,  or  other  data  is  not  to  be  regarded 
by  implication  or  otherwise  as  in  any  manner  licensing  the  holder  or  any  other 
person  or  corporation  or  conveying  any  rights  or  permission  to  manufacture,  use 
or  sell  any  patented  invention  that  may  in  any  way  be  related  thereto. 

This  report  has  been  reviewed  by  the  Information  Office  (01)  and  is 
releasable  to  the  National  Technical  Information  Service  (NTIS).  At  NTIS, 
it  will  be  available  to  the  general  public,  including  foreign  nations. 

This  technical  report  has  been  reviewed  and  is  approved  for  publication. 


FOR  THE  COMMANDER 


Technical  Director,  ARTO 
Project  Officer 


z /\  / -s  yy 

‘vuONALD  L.  LAMBERSON 


Colonel,  USAF 

Chief,  Advanced  Radiation  Technology 
Office 


/ / *"  -i 

Vv 

V 


DO  NOT  RETURN  THIS  COPY.  RETAIN  0P.  DESTROY. 


SECURITY  Cl  ASSl  PiCATiQn  Of  This  page  'Whan  Data  Entarad ) 


EPORT  DOCUMENTATION  PAGE 


2-  GOVT  ACCESSION  NO 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3 »ECi^iEnt  s Catalog  number 


* TITLE  rrnnd  Subtnl*) 


UNSTABLE  RESONATOR ^THEORY* 


-•  fl^cully,  S.  J^^oore/^FL  J^icCarthy 

— ^ 
Optical  Sciences  Center 
University  of  Arizona 

Tucson,  Arizona  85721 

M.  CONTROLLING  OFFICE  NAME  AnO  AOORE5S 

Air  Force  Weapons  Laboratory  (AR)  ( 

Kirtland  Air  Force  Base,  NM  87117 


Ijjjj 


3 CONTRACT  QP  GRANT  NuM3£Rf»; 

| F296/h-75-C-O08 8/ 


10.  PROGRAM  ELEMENT.  PROJECT,  task 
AREA  & *OR<  UNIT  NUMBERS 


6261 00/31 7J01S9 


Mar#  077 


^4.  MONITORING  AGENCY  NAME  A AOORESSfif  dlffaeant  Irom  Controlling  Office;  1 IS.  SECURITY  CLASS,  (ol  this  r«portj 


(p)  '/Ip 


UNCLASSIFIED 


I Sr.  OECLASSlPlCATlON  DOWNGRADING 
SChEOULS 


DISTRIBUTION  STATEMENT  'of  this  Raport, 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (ol  tha  abstract  ant  a rad  in  Block  20,  it  dltlarant  from  Raport) 


\ it.  supplement  ary  notes 


1 19.  KEY  WORDS  (Contlnua  on  raaaraa  aida  it  nacaasary  and  idantlty  by  block  numbar ) 


Unstable  Resonators 
Laser  Theory 
D1 ffraction 


^ABSTRACT  rContlnua  on  raaaraa  aida  It  nacaaaary  and  ioantlty  by  block  numbar ) 


This  report  on  research  on  unstable  resonator  theory  consists  of  three  parts. 
Section  I gives  an  overview  of  the  work.  In  Section  II  a simple  model  is 
described  for  estimating  the  intensity  distributions  inside  and  emerging  from 
lasers  with  high-Fresnel-number  unstable  resonators.  In  Section  III  modes  of 
strip  confocal  unstable  resonators  containing  spatially  dependent  gain  are 
determined  by  means  of  a theory  that  is  practical  even  at  large  Fresnel 
numbers. 


DO  1 janM73  1473  \eO.TIOn  op  I NOV  «s  is  OBSOLETE 


UNCLASSIFIED 


SECURITY  CL  ASSI  PiC  ATion  op  This  PAGE  (Whan  Data  Entarad) 


& 


CONTENTS 


Section 

I 

II 

III 


Page 

WORK  OVERVIEW  1 

LASERS  WITH  UNSTABLE  RESONATORS  IN  THE  GEOMETRIC 
OPTICS  LIMIT  3 

THEORY  OF  MODES  IN  A LOADED  STRIP  CONFOCAL 

UNSTABLE  RESONATOR  34 

REFERENCES  87 


ILLUSTRATIONS 


Figure  Page 

Fig.  1.  Geometry  of  the  symmetric  unstable  resonator.  17 

Fig.  2.  Interior  and  exit  rays  in  the  symmetric  unstable  resonator.  18 

Fig.  3.  Qualitative  transverse  dependence  of  intensity  in  the  19 

output-mirror  plane  in  the  geometrical  limit. 

Fig.  4.  Symmetric  resonator.  M = 2.5.  GqD  = 4 for  0 < z < 0. 

a)  On-axis  symmetry,  b)  Output-mirror-plane  Intensity.  20-21 

Fig.  5.  Symmetric  resonator.  M 3 2.5.  GqD  = 4 for  D/4  < z < 3D/4. 

On-axis  intensity.  22 

Fig.  6.  Geometry  of  the  confocal  unstable  resonator.  23 

Fig.  7.  Confocal  resonator.  M * 2.5.  G D 3 4 for  Q < z < d. 

3 o 

a)  On-axis  intensity,  b)  Output-mirror-plane  intensity.  24-25 

Fig.  8.  Confocal  resonator.  M * 2.5.  GQ0  3 4 for  z < D/2. 

a)  On-axis  intensity,  b)  Output-mirror-plane  intensity.  26-27 

Fig.  9.  Confocal  resonator.  M * 2.5.  GQD  3 4 for  D/4  < z < 3D/4. 

On-axis  intensity.  28 

Fig.  10.  Confocal  resonator.  M 3 1.5.  GQD  3 4 for  0 < z < D. 

On-axis  intensity.  29 

Fig.  11.  Noniform  grid  used  in  calculations  where  GQ  depends  on  r. 

The  symbol  X on  the  horizontal  axis  marks  the  focal  point.  30 

iii 


Figure  Page 

Fig.  12.  Confocal  resonator.  GoD  = 4 exp  [-0.5(r/a  - 0 . 75 ) 2 ] for 
r > 0.75a.  Output-mirror-plane  intensity. 

a)  M = 1.5  b)  M = 2.5.  31-32 

Fig.  13.  Geometry  of  the  confocal  unstable  resonator.  71 

Fig.  14.  Geometry  of  the  strip  confocal  unstable  resonator  is  illustrated 
for  a resonator  with  M = 1.5.  The  x coordinate  is  shown  on  a 
much  expanded  scale  relative  to  the  o coordinate.  Edge  dif- 
fraction effects  are  accounted  for  in  terms  of  cylindrical  waves 
emanating  from  points  (pn,+xn),  where  pn  = dir  and  xn  = aMn. 

The  points  (p0>xo)  and  (p^  ,x-| ) act  as  edges  of  apertures. 

Some  of  the  shadow  boundaries  on  one  side  of  the  resonator 

are  shown  with  dotted  lines.  For  a resonator  with  large  M,  the 

shadow  boundaries  would  not  come  so  close  to  the  mirror  edges.  72 

Fig.  15.  Ray  paths  along  which  the  gain  is  integrated  to  give  Tn(+a) 
are  shown  for  n = 1,2, 3, 4.  The  output  mirror  edge  is  in  the 
upper  right  corner  of  each  diagram  and  the  resonator  axis  is 
along  the  bottom.  Rays  below  the  axis  (not  shown)  are  the 
mirror  images  of  those  above.  Those  rays  which  return  to  x = a 
do  so  by  reversing  the  path  by  which  they  zig-zag  toward  the 
axis.  The  diagrams  are  shown  for  a resonator  with  M = 2.5.  73 

Fig.  16.  Coordinate  systems  used  to  evaluate  the  integral  I in  Eg.  (50).  74 

Fig.  17.  The  functions  Tn'(x)  for  n * 1 through  7 are  shown  for  a resonator 
with  M = 2.9  containing  a Gaussian  gain  distribution.  The  upper 
curve  is  for  n = 1.  As  n increases  Tn'(x)  approaches  a symmetric 
limit  function.  75 

iv 


i 


Figure  Page 

Fiq.  18.  Modulus  of  A vs.  effective  Fresnel  number  for  M = 2.9,  symmetric 
roots,  a)  Bare  resonator,  b)  Resonator  containing  Gaussian  gain 
distribution.  76-77 

Fiq.  19.  Modulus  of  A vs.  effective  Fresnel  number  for  M = 2.0,  symmetric 
roots,  a)  Bare  resonator,  b)  Resonator  containing  Gaussian  gain 
distribution,  c)  Resonator  containing  negative  Gaussian  gain  distri- 
bution. 78-30 

Fiq.  20.  Phase  and  intensity  of  the  lowest-loss  mode  of  a bare  resonator, 

for  Feff  * 16.4,  M = 2.9,  u = 0.96145  - 0.0042704i.  81 

Fig.  21.  Phase  and  intensity  of  the  next-to-lowest-loss  symmetric  mode 
of  a base  resonator,  for  Fgff  = 16.4,  M » 2.9,  y = 1.1001 
+ 1 . 1 502i . 82 

Fig.  22.  Phase  and  intensity  of  the  lowest-loss  mode  of  a resonator 

containing  a Gaussian  gain  distribution,  for  F ,,  = 16.4, 

ett 

M = 2.9,  y * 0.97999  - 0.0028037i.  83 

Fig.  23-  Phase  and  intensity  of  the  lowest-loss  mode  of  a bare  resonator, 

for  Ffiff  = 8.4,  M = 2.0,  y = 0.95172  - 0.00602461.  84 

Fig.  24.  Phase  and  Intensity  of  the  lowest-loss  mode  of  a resonator 

containing  a Gaussian  gain  distribution,  for  Feff  * 8.4,  M = 2.0, 

y « 0.9716  - 0.00388151.  85 

Fig.  25.  Phase  and  intensity  of  the  lowest-loss  mode  of  a resonator 
containing  a negative  Gaussian  gain  distribution,  for  Fg^  * 

8.4,  M * 2.0,  y » 0.93869  - 0.00721241.  86 


v/vi 


SECTION  I 


WORK  OVERVIEW 

This  report  summarizes  work  completed  on  unstable  resonator  theory  for  the 
BOA  contract  F29601-75-C-0088  between  the  Air  Force  Weapons  Laboratory,  Kirtland 
Air  Force  Base,  and  the  University  of  Arizona  Optical  Sciences  Center.  The  first 
section  serves  as  an  overview  of  the  work.  The  second  and  third  sections  present 
preprints  of  two  papers  on  unstable  resonator  theory  to  be  published,  probably  in 
the  Journal  of  the  Optical  Society  of  America. 

1.  Geometrical  Theory  of  a Loaded  Resonator 

We  have  developed  a simple  geometrical  model  for  estimating  the  intensities 
inside  of  and  emerging  from  high-Fresnel-number  unstable  resonators.  This  model 
uses  first-order  rate  equations  to  determine  the  wave  amplitudes.  The  wavefronts 
are  taken  to  be  the  same  as  in  Siegman's  early  geometrical  theory.  The  amplitudes 
in  the  core  of  the  resonator  (oscillator  region)  are  found  by  solving  a first-order 
differential  equation  along  the  resonator  axis  subject  to  boundary  conditions  on 
the  mirrors.  The  amplitudes  in  the  outer  part  of  the  resonator  (amplifier  region) 
are  found  by  solving  the  rate  equations  in  the  outer  region  and  matching  their 
solutions  to  the  core  solution.  The  field  is  assumed  to  be  chopped  off  discon- 
tinuously  at  the  shadow  boundaries.  It  the  unsaturated  gain  GQ  depends  on  the 
transverse  coordinate,  then  a simple  iterative  procedure  is  used  to  determine  the 
self-consistent  steady-state  solution.  Although  the  predictions  of  the  geometrical 
theory  are  crude,  the  theory  is  useful  because  it  is  simple  and  requires  very  little 
computer  time  to  perform  calculations.  This  theory  is  described  in  the  preprint 
of  Section  II,  "Lasers  with  Unstable  Resonators  in  the  Geometrical  Optics  Limit." 
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2.  Theory  Including  Edge  Diffraction  as  well  as  Transverse  Gain  Distribution 

We  have  developed  a theory  that  enables  one  to  use  only  a small  amount  of 
computer  time  to  calculate  modes  of  loaded  unstable  resonators  with  sharp-edged 
mirrors.  The  theory  has  been  developed  so  far  only  for  strip  (two-dimensional) 
resonators,  but  it  appears  that  it  can  be  generalized  to  resonators  with  rectangu- 
lar mirrors  in  a straightforward  manner.  Besides  the  advantage  of  ease  of  compu- 
tation, the  theory  offers  other  advantages  over  the  fast  Fourier  transform  (FFT) 
method.  It  is  applicable  at  high  effective  Fresnel  numbers.  Indeed,  it  is 
expected  to  become  more  accurate  as  Fgff  increases.  Calculation  for  F ^ of  several 
hundred  is  feasible.  On  the  other  hand,  it  also  works  well  for  low  Fef<r  (near 
unity)  provided  that  the  magnification  M is  not  too  small.  The  theory  can  be  used 
to  calculate  any  of  the  resonator  modes,  not  just  the  lowest-loss  mode.  However, 
the  theory  does  not  include  effects  of  inhomogeneities  in  the  index  of  refraction. 
Also,  it  is  not  clear  whether  an  analogous  theory  can  be  developed  for  resonators 
with  round  mirrors. 

This  theory  contains  elements  of  P.  Horwitz's  "Asymptotic  Theory  of  Unstable 

Resonator  Modes, elements  of  our  geometrical  theory,  and  elements  of  Rubinowicz's 

? 

line-integral  approach  to  edge  diffraction.  As  in  Horwitz's  theory,  the  self- 
consistency  condition  on  the  gain  is  found  by  solving  for  the  roots  of  a polynomial. 

This  theory  is  described  in  the  preprint  of  Section  II,  "Theory  of  Modes  in 
a Loaded  Strip  Confocal  Unstable  Resonator."  A talk  on  this  theory  was  given  at 
Kirtland  on  April  26,  1976. 
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SECTION  II 


LASERS  WITH  UNSTABLE  RESONATORS  IN  THE  GEOMETRICAL  OPTICS  LIMIT 
1.  Introduction 

To  develop  some  physical  understanding  of  the  behavior  of  lasers  with  unstable 
resonators  and  to  obtain  numerical  estimates  of  intensities,  it  is  useful  to  con- 
sider models  which  are  relatively  simple.  This  treatment  uses  geometrical  optics 
and  rate  equations  to  describe  the  single-frequency  operation  of  a laser  containing 
a medium  with  saturable  gain.  The  treatment  is  valid  for  resonators  with  high 
effective  Fresnel  number.  However,  since  diffraction  of  light  at  the  mirror 

edges  is  neglected,  the  high-frequency  low-amplitude  ripples  in  the  intensity 

1 -4 

which  are  found  in  more  exact  treatments  do  not  appear  in  the  present  theory. 
Also,  the  intensity  drops  to  zero  discontinuously  at  the  shadow  boundaries.  These 
defects  are  offset  by  the  ease  with  which  calculations  can  be  made  using  our 
theory.  The  more  exact  theories  require  time-consuming  and  expensive  computation. 
Our  theory  allows  a rapid  determination  of  the  resonator  parameters  which  are  likely 
to  give  optimal  results  in  a more  exact  calculation. 

In  the  simple  case  where  the  unsaturated  gain  GQ  is  independent  of  the  trans- 
verse coordinates  our  model  allows  the  calculation  of  the  field  distribution 
inside  the  laser  by  solving  a first-order  ordinary  differential  equation  for  the 
behavior  in  the  core  region  and  matching  this  to  the  solution  of  a similar  equation 
describing  the  "exit  rays"  leaving  the  laser.  If  GQ  depends  on  the  transverse 
coordinates,  then  a simple  iterative  calculation  can  be  used  to  determine  the  off- 
axis  field.  Our  model  is  sufficiently  general  to  describe  situations  where  the 
dynamics  of  the  active  medium  breaks  the  axial  symmetry  of  the  laser  (for  instance, 
a flowing  medium).  Indeed,  the  optical  rate  equations  restrict  the  present 
analysis  to  the  case  where  the  laser  and  the  field  which  it  contains  are  axially 
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symmetric.  We  also  neglect  all  refraction  and  diffraction  and  assume  that  the 
resonator  is  perfectly  aligned.  We  discuss  two  geometries,  the  symmetric  resona- 
tor and  the  confocal  unstable  resonator.  The  initial  discussion  is  in  terms  of 
the  symmetric  resonator.  Numerical  solutions  for  the  intensities  are  presented 
for  both  geometries,  but  with  most  attention  being  given  to  the  confocal  unstable 

resonator.  The  geometrical  aspects  of  our  treatment  are  based  on  early  work  of 
5 

Siegman.  The  rate  equations  which  we  use  are  the  counterpart  for  the  unstable 
resonator  of  equations  used  by  Rigrod^  to  analyze  one-dimensional  lasers. 

2.  Basic  Rate  Equations 

5 

Siegman,  in  his  first  paper  on  unstable  resonators  , considered  that  the 
field  inside  an  unstable  resonator  consists  of  two  spherical  waves  which,  for  the 
symmetric  resonator  pictured  in  Fig.  1,  emanate  from  virtual  centers  C-|  and  C2- 
The  positions  of  these  centers  are  determined  by  the  resonator  length  D and  the 
mirror  curvatures.  It  is  easy  to  show  that  this  pair  of  waves  reflect  into  each 
other  if  the  mirrors  have  the  form  of  hyperboloids  with  foci  at  C-j  and  C2>  Thus, 
for  the  symmetric  resonator  we  take  the  equations  of  the  left  and  right  mirrors 
to  be  p2  = p-|  + D and  = p2  + D,  where  p^ , and  p2,  e2  are  spherical  coordi- 
nates centered  at  and  C2-  We  make  a paraxial  approximation  where  necessary  by 
replacing  sin  e..  by  and  cos  e..  by  1 . In  this  approximation  it  does  not  matter 
whether  the  mirrors  are  hyperboloidal , spherical,  or  whatever,  as  long  as  the 
mirror  curvatures  remain  the  same.  The  paraxial  approximation  is  an  excellent 
one  in  many  cases  of  practical  interest,  but  could  be  eliminated  at  the  cost  of 
increased  mathematical  complexity.  If  both  mirrors  of  the  symmetric  resonator 
have  radius  a and  L is  the  distance  between  C-|  and  C2,  we  require  that  the  outer 
angle  2a/(L-D)  of  the  cones  of  radiation  emerging  from  the  resonator  be  much  less 
than  one. 
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It  is  easy  to  see  that  if  one  combines  two  spherical  waves  which  are  solu- 
tions of  the  free-space  Helmholtz  equation,  then  one  can  match  the  phases  of 
the  waves  on  the  mirrors  by  taking  the  wave  vector  k = mr/D.  However,  the  ampli- 
tudes do  not  match.  This  is  not  surprising,  since  we  do  not  expect  a steady 
state  in  this  case.  However,  we  can  easily  match  the  amplitudes  if  we  introduce 
a medium  with  a linear  gain  y into  the  resonator  and  choose  y so  as  to  just 
offset  the  geometrical  losses.  Then  the  field  becomes 

E “ p-|  1 exp  (imrp^/D  + yp-j)  + (-l)n  1 exp  (int^/D  + YP2).  (1) 

This  is  zero  on  the  mirrors  within  the  paraxial  approximation  if  we  take 
y = (1/D)  1 M vhere  M = (L+D)/(L-D)  is  the  one-way  magnification.  The  waves 

first  decrease  in  amplitude  because  of  geometrical  spreading  as  one  moves  away 
from  their  centers  and  then  increase  in  amplitude  again  because  of  the  gain.  In 
this  way  the  amplitudes  can  be  made  to  match  on  the  mirrors. 

We  expect  that  something  of  this  sort  also  happens  with  a nonlinear  medium 
in  the  resonator.  However,  for  the  case  of  a nonlinear  medium  above  threshold, 
the  intensity  adjusts  itself  so  that  the  gain  G equals  the  loss  in  the  steady 
state.  We  anticipate  that  G will  be  a function  of  position  and  that  the  amplitude 
of  the  spherical  waves  may  depend  on  e..  as  well  as  p • . 

To  obtain  the  basic  rate  equations  of  our  theory,  we  begin  with  the 
wave  equation 

V2E  = c'232E/3t2  - 2c_1G3E/3t  . (2) 

The  saturated  gain  G has  been  put  into  the  wave  equation  in  the  form  of 
a negative  conductivity.  The  gain  is  a function  of  position,  but  we 
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assume  that  the  index  of  refraction  is  independent  of  position.  This 
can  be  a valid  assumption  if  the  laser  operates  in  the  center  of  the 
gain  profile.  We  take  each  spherical  wave  inside  the  laser  to  have  the 
form 

E = [f(p,e,t)/p]  eikp"ioyt  , (3) 


where  k » G.  The  function  f is  assumed  to  be  slowly  varying  in  both 
space  and  time.  If  we  insert  Eq.  (3)  into  Eq.  (2)  and  put  oi  = ck,  we 
obtain  the  following  equation  for  f: 


!i+  2ik  11  + -LJ- 

2 3o  2-  36 


3f  \ 


:-2fl 


2ioj 


'•at 


3f 

at 


- 2c'1G(§£  - iuf) 


(4) 


2 2 2 2 
If  we  neglect  3 f/ap  compared  to  2ik3f/3p,  neglect  3 f/3t  compared  to 

2ia)3f/3t,  and  neglect  3f/3t  compared  to  iwf,  Eq.  (4)  reduces  to  the  following 

paraxial  wave  equation: 


2ik 


‘g ♦ 


3f 

at 


- Gf)  + 


p2e 


(5) 


Equation  (5)  is  sufficiently  general  to  handle  edge-diffraction  effects. 
However,  we  wish  to  make  a further  approximation  to  Eq.  (5)  by  dropping 
the  term  involving  derivatives  with  respect  to  9.  This  gives  the  rate 
equation 


3f 

9p 


+ 


Gf 


(6) 


Equation  (6)  can  be  used  to  describe  time-varying  phenomena  such  as  the 


6 


build-up  of  the  laser  field  when  the  pumping  is  switched  on.  In  the  steady  state 
3f/3t  = 0,  so  that 


3f/3P  = Gf  . (7) 

Equation  (7)  is  the  basic  rate  equation  of  our  theory. 

The  neglect  of  the  term  in  Eq.  (5)  containing  derivatives  with  respect  to  9 

represents  a great  simplification  of  the  problem,  but  eliminates  the  possibility 

of  accounting  for  edge-diffraction  effects.  The  rays  associated  with  Eq.  (7)  all 

pass  through  the  point  p * 0,  and  there  is  no  possibility  of  ray  crossing.  Edge 

diffraction  gives  rise  to  rays  which  emanate  from  the  mirror  edges  and  bounce 

back  and  forth  across  the  resonator,  crossing  in  a complicated  way.  The  wavelength 

dependence  of  the  resulting  diffraction  pattern  is  conveniently  parameterized  in 

terms  of  the  effective  Fresnel  number  F ^ (also  sometimes  called  the  equivalent 

Fresnel  number  N ).  One  might  think  that  for  large  effective  Fresnel  number  Fg^ 

(i.e.,  large  k)  this  approximation  would  always  be  a good  one.  Actually  the 

situation  is  more  complicated.  As  k increases,  the  transverse  oscillations  in  f 

become  more  and  more  rapid,  and  it  is  not  immediately  clear  that  the  neglected 

term  is  small.  Numerical  calculations  of  Horwitz1  show  that  several  modes  can 

exist  within  the  resonator  and  that  only  the  lowest-loss  mode  approaches  the 

result  expected  from  geometrical  theory  as  k -*•  » . This  mode  does  not  approach 

the  geometrical  limit  in  a smooth  manner,  but  rather  it  develops  more  and  more 

rapid  but  smaller  and  smaller  amplitude  ripples  as  k increases.  The  phase  of  f 

is  nearly  constant  in  the  resonator  core  even  at  fairly  low  F but  the  ripples 

in  the  modulus  of  f die  out  rather  slowly  as  F ^ increases.  They  are  about  :15* 

of  the  mean  value  for  F^  near  50.  The  magnitude  of  the  ripples  can  be  made  less 
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than  this  by  using  mirrors  with  tapered  reflectivity  profiles.  ’ The  calculations 
of  Horwitz  apply  to  resonators  with  rectangular-edged  mirrors  or  to  strip 
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resonators  with  round 


resonators.  Unfortunately  no  calculations  for  high  F ^ 
mirrors  have  so  far  been  carried  out.  Probably  the  main  difference  to  be  expected 
with  round  mirrors  is  an  increase  in  intensity  in  a region  very  close  to  the 
resonator  axis.  This  increase  results  from  constructive  interference  of  diffracted 
light  coming  from  the  mirror  edges.  The  size  of  the  region  should  shrink  toward 
zero  as  Fg^  increases.  In  summary,  it  appears  that  Eq.  (7)  is  valid  in  an 
average  sense  for  the  lowest-loss  mode  at  high  Fg^. 

Figure  2 shows  ray  paths  for  the  symmetric  resonator.  We  see  that  the  reso- 
nator is  divided  into  two  regions.  In  the  core  the  left-  and  right-travel ing 
rays  are  both  present  and  both  contribute  toward  saturating  the  gain  medium.  In 
the  outer  part  of  the  resonator  the  "exit  rays"  leaving  the  resonator  do  not  over- 
lap. Let  us  consider  the  change  in  intensity  along  the  exit  rays.  For  a homo- 
geneously broadened  medium  we  may  take  G to  be  a function  of  the  intensity  I 
according  to 


G - G0(l  + I/I0)_1  , (8) 

where  GQ  is  the  unsaturated  gain,  which  may  be  taken  to  be  constant  or  a given 
function  of  position,  and  I is  the  saturation  intensity.  We  choose  f real  and 

9 9 

normalize  the  field  so  that  I = f /p  . If  we  insert  Eq.  (8)  into  Eq.  (7),  we  get 

3f/3p  = G0f/[1  + f2/(I0P2)]  • (9) 

This  equation  must  be  solved  numerically,  although  asymptotic  properties 
of  the  solution  are  easily  derived.  For  instance,  if  Gq  is  constant  and 
p is  large,  then 

f * (2GqI0/3)1/2  p3/2  , (10) 
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so  that  the  intensity  increases  linearly  with  p for  large  p. 


3.  Symmetric  Resonator  Intensities  for  the  Case  That  S.  Does  Mot 
Depend  on  r u 

We  write  the  field  in  the  symmetric  resonator  in  the  form 

E 3 of1  eikpi  f(p1>9l)  + (-Dn  p2_1  eikp2  g(p2,e2),  (11) 

which  is  a generalization  of  Eq.  (1).  The  functions  f and  g are  assumed  to  be 
slowly  varying.  In  order  that  E be  zero  on  the  mirrors,  we  require  that 

f(p1.01)/p1  3 g(p2,e2)/p2  (12) 

on  the  mirrors.  Using  the  geometrical  relationship  between  the  p^  ,9^  and  p2,e2 
coordinates  on  the  mirrors,  Eq.  (12)  becomes 


f[**(L+D),97]  = 

Mg[>s(L-D)  ,M97] , 

(13) 

Mfft(L-D),M92] 

= g[)s(L+D),92]. 

(14) 

These  relations  could  be  easily  generalized  if  the  mirrors  were  not  perfectly 
reflecting,  but  had  reflectivities  which  were  functions  of  the  radial  coordinate 
on  the  mirrors. 

Let  us  suppose  that  in  some  region  suitably  close  to  the  axis  it  is  a good 
approximation  to  assume  that  GQ  is  independent  of  r.  Then  in  this  region  it  is 
consistent  to  assume  that  the  field  itself  is  independent  of  r,  so  that  f does 
not  depend  on  9-j  and  g does  not  depend  on  92>  Using  Eq.  (7),  we  have 

df/dp^  3 
dg/dp2  = 
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Gf  , 
Gg  . 


(15) 

(16) 


In  the  paraxial  approximation  p2  = L - p^ • We  take  the  Wronskian  of  Eqs.  (15) 
and  (16)  by  multiplying  (15)  by  g and  (16)  by  f and  subtracting.  This  gives,  upon 
integrating. 


f(p^)  g(p2)  = C , 


(17) 


where  C is  a constant. 

Using  the  Cartesian  z coordinate  shown  in  Fig.  1,  we  define  fQ  and  H(z)  by 
the  equations 

f£*s(L-D)  + z]  = f eH(z)  , (18) 

u < 

H(o)  = 0 . (19) 

It  follows  from  Eqs.  (15)  and  (18)  that 

dH/dz  = G = G0/(I  + I/I0)  . (20) 

The  total  time-averaged  intensity  I is  found  from  Eq.  (10)  to  be 

1 = P]  2 f2(p-|)  + p2  2 92(pj)  ‘ (*1)n  2P-|-1  P2"1  9(p2) 

x cos  [nit  (P]  - p2)/D]  . (21) 


The  last  term  is  a rapidly  oscillating  function  of  position  which  averages  to 
zero,  and  it  will  be  neglected.  We  wish  to  write  Eq.  (20)  as  a differential 
equation  for  H(z).  Combining  Eqs.  (17),  (13),  (14),  and  (18),  we  have,  at  the 
left  and  right  mirrors. 


f2  [»s(L-D)]  M ■ fQ2  M , 

(22) 

f2  [«*(L+D)]/M  = fQ2  e2H(D)/M  . 

(23) 
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By  combining  Eqs.  (22)  and  (23),  we  obtain  the  boundary  condition 

H(D)  = IH  M . (24) 

Also,  using  Eqs.  (18)  and  (22)  in  Eq.  (17),  we  have 

g[»s(L+D)  - z]  - M fQ  e'H(z)  . (25) 

Using  Eqs.  (21,  (18),  and  (25)  in  Eq.  (20),  we  obtain  the  desired  differential 
equation  for  H(z) : 

dH/dz  = Gq/{  1 + [(>s(L-D)  + z)'2  e2H  + M2(^(L+D)  - z)-2  e‘2H]  fQ2/I0>  (26) 

This  equation  must  be  solved  subject  to  the  boundary  conditions  (19)  and  (24). 

2 

Both  of  these  boundary  conditions  can  be  satisfied  only  if  f /I  has  a particular 
value.  This  value  determines  the  steady-state  intensity.  Solutions  exist  only 
if  the  laser  is  above  threshold  (G0D  > m M). 

If  Gq  is  independent  of  r for  all  r,  then  the  intensities  found  by  solving 
Eq.  (26)  are  valid  out  to  the  region  of  the  exit  rays.  The  intensities  along  the 
exit  rays  may  be  found  by  matching  solutions  of  Eq.  (9)  to  the  interior  solution 
at  the  conical  boundaries  9.  = 2a/ (L-D)  separating  the  interior  from  the  exit  ray 
region  (see  Fig.  2).  Since  the  medium  is  less  saturated  in  the  exit  ray  region 
than  in  the  interior,  we  expect  the  intensity  to  increase  with  angle  within  the 
exit  cones  out  to  the  "cut-off"  angle  3 2a(L-D),  where  I drops  to  zero,  as 
shown  qualitatively  in  Figure  3. 

In  Figures  4 through  6 we  show  numerical  solutions  for  the  symmetric  unstable 
resonator  in  which  GQ  is  independent  of  r.  Here  we  assume  that  M = 2.5  and  that 
the  mirrors  have  equal  radii.  Figures  4a,  5,  and  6 show  the  intensities  1^  and 
I2  (in  the  core  of  the  resonator)  of  the  waves  traveling  to  the  right  and  left. 
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respectively.  These  are  given  in  units  of  I and  plotted  as  a function  of  z/D. 
Computed  points  associated  with  I-j  and  are  indicated  by  +'s  and  x's  respec- 
tively. Figure  4b  shows  the  intensity  on  either  mirror  plane  as  a function  of 
r/z  in  the  domain  1 < r/a  < M.  For  r/a  < 1 the  intensity  is  independent  of  r. 
For  r/a  > M the  intensity  drops  discontinuously  to  zero.  Figure  4 is  for  the 
case  where  GQD  = 4 throughout  the  cavity.  Figure  5 is  for  the  case  where 
GqD  = 4 for  D/4 < z < 3D/4  and  zero  otherwise. 

4.  Confocal  Unstable  Resonator  Intensities  for  the  Case  That 
G«  Does  Not  Depend  on  r 

The  confocal  unstable  resonator  is  the  configuration  of  choice  in  most 
applications  because  the  output  emerges  in  the  form  of  a collimated  beam.  The 
geometry  of  the  confocal  unstable  resonator  is  shown  in  Figure  6.  The  field 
consists  of  a plane  wave  moving  to  the  right  and  a spherical  wave  moving  to  the 
left  and  centered  on  a point  C a distance  d to  the  right  of  the  convex  mirror. 
The  cavity  length  is  D.  The  magnification  of  the  spherical  wave  in  traveling 
from  the  right  mirror  to  the  left  mirror  is  M = (D+d)/d.  The  left  and  right 
mirrors  may  be  taken  to  be  parabaloids  given  by  the  equations  p = z + D + d and 
p = z - D + d.  The  field  in  the  cavity  has  the  form 

E = eik(z+D+d)  f (z,r)  - p-1  eikp  g(p,e)  , (27) 

where  k = nir/D.  In  order  that  E be  zero  on  the  mirrors,  we  require  that 

f(0,r)  = (D+d)'1  g[D+d,r/(D+d)]  , (28) 

f (D,r)  = d"1  g(d,r/d)  . 

The  amplitudes  f and  g obey  the  rate  equations 


3f(z,r)/3Z 

= Gf(z.r)  , 

(30) 

ag(p,e)/3p 

= Gg(p,e)  . 

(31) 
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We  assume  that  the  concave  mirror  is  large  enough  so  that  all  of  the  radiation 
comes  out  to  the  right.  If  GQ  is  essentially  independent  of  r in  a core  region 
suitably  close  to  the  axis,  then  the  field  may  be  taken  to  be  independent  of  r in 
this  core  region.  If  the  core  region  extends  out  to  the  cone  0 = a/d,  which  is 


the  boundary  between  the  interior  of  the  laser  and  the  exit-ray  region,  then  the 
solution  in  the  interior  may  be  extended  out  to  the  shadow  boundary  r = Ma  by 
means  of  the  exit-ray  equation 

3f/3z  = GQf/[l  + f2/Io]  . 

(32) 

This  equation  is  found  by  using  Eqs.  (30)  and  (8). 

The  differential  equation  giving  the  intensities  in  the  core  region  may  be 
derived  by  the  same  method  that  was  used  for  the  symmetric  resonator.  Taking 
the  Wronskian  of  Eqs.  (30)  and  (31),  we  find  that 

f(z)  g(p)  = C , 

(33) 

where  C is  a 

constant.  We  define  H(z)  and  f by  the  equations 

f(z)  = fQ  eH(z)  , 

(34) 

H(0)  = 0 . 

(35) 

Combining  Eqs 

. (33),  (28),  (29),  and  (34),  we  have,  at  the  left 

and  right  mirrors 

C = (D+d)  f2(0)  = (D+d)  fQ2  , 

(36) 

C = d f2( D)  = d fQ2  e2H(D)  . 

(37) 

By  combining  Eqs.  (36)  and  (37),  we  obtain  the  boundary  condition 

H(D)  = H in  M . (38) 
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Note  that  the  intensity  at  the  right  end  of  the  cavity  is  a factor  of  M greater 
than  the  intensity  at  the  left  end.  Using  Eqs.  (34)  and  (36)  in  Eq.  (33),  we 
have 


g(D  + d - z)  = (D+d)  fQ  e~H(z)  . (39) 

Using  Eqs.  (34),  (39),  and  (8),  Eq.  (30)  may  be  written  as  a differential  equation 
for  H(z) : 

dH/dz  = g0/{l+[e2H  + [ ( D+d ) / ( D+d-z) ] 2 e'2H]  fQ2/Io>  . (40) 

This  equation  must  befsolved  subject  to  the  boundary  conditions  (35)  and  (38). 

In  Figures  7 through  10  we  show  numerical  solutions  for  the  confocal  unstable 
resonator  in  which  GQ  is  independent  of  r.  Figure  7 is  for  a resonator  with 
M = 2.5  and  GQD  = 4.  Figure  7a  shows  the  intensities  in  the  core,  while  Figure 
7b  shows  the  intensity  in  the  plane  of  the  convex  mirror.  Figures  8 and  9 are 
also  for  a resonator  with  M = 2.5  and  GQD  = 4,  except  that  the  gain  medium  occupies 
only  the  regions  0 < z<  0/2  and  D/4  < z < 3D/4,  respectively.  Figure  10  is  for 
a resonator  with  GQD  = 4 everywhere,  but  with  M = 1.5. 

5.  Confocal  Unstable  Resonator  Intensities  for  the  Case  That  G Does 
Depend  on  r 0 

In  general  it  is  not  a good  approximation  to  assume  that  GQ  is  independent  of 
r all  the  way  out  to  the  shadow  boundaries.  In  this  section  we  show  how  intensi- 
ties in  the  confocal  unstable  resonator  may  be  calculated  for  an  arbitrary  spatial 
dependence  of  Gq.  For  definiteness  we  assume  that  the  core  region  inside  of 
which  Gq  may  be  considered  independent  of  r is  r < a/M.  Smaller  core  sizes 
could  be  used  without  much  extra  work.  Inside  the  core  the  solutions  found  in 
Section  IV  are  still  valid.  The  solutions  outside  the  core  may  be  found  by 
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integrating  che  rate  equations  (30)  and  (31)  along  the  light  rays  as  they  bounce 
out  of  the  resonator,  matching  the  intensities  of  the  left-  and  right-travel ing 
waves  on  the  mirror  planes.  However,  the  saturated  gain  is  a function  of  the 
intensities  1^  and  ^ which  are  not  known  until  the  calculation  is  completed. 
Moreover,  the  use  of  the  Wronskian  to  cast  the  problem  into  the  form  of  a 
differential  equation  is  not  possible  if  f depends  on  r and  g depends  on  e.  We 
overcome  this  difficulty  by  doing  an  iterative  calculation  in  which  the  intensi- 
ties and  the  gain  are  calculated  successively  until  both  converge  to  a self- 
consistent  solution.  We  find  that  this  convergence  occurs  after  about  five 
iterations  and  takes  very  little  computer  time. 

In  order  to  carry  out  the  numerical  integration,  I-j,  Ig,  GQ,  and  G must  be 

evaluated  on  a two-dimensional  grid.  We  use  the  non-uniform  grid  (z  ,r  ),  where 

j 3 v m n 

zm  = [DM/(M-1)]  (1-M'm/N),  n=0,...,N,  (41) 

rn  = aMn/N_1  , n = 0.....2N  . (42) 

This  grid  is  illustrated  in  Figure  11  for  N = 10.  The  advantage  of  this  grid 
is  that  both  the  right-  and  left-traveling  rays  pass  through  the  grid  points,  so 
that  no  interpolation  is  necessary.  We  integrate  using  the  trapezoidal  rule, 
stepping  up  one  value  of  n at  a time  from  the  known  solution  in  the  core,  first 
calculating  l£  and  then  I,.  It  is  possible  to  do  the  iteration  at  each  value 
of  n before  proceeding  to  the  next  value.  This  would  be  convenient  if  the  calcu- 
lation were  to  be  performed  on  a desk  calculator  with  limited  memory. 

Figure  12  shows  some  results  of  calculations  using  our  iterative  scheme. 

We  take  GQ  here  to  depend  on  r according  to 

Gq0  = 4 exp  [-0. 5(r/a  - 0.75)2]  (43) 

for  r > 0.75a.  Figures  12a  and  12b  show  the  intensity  on  the  right  mirror  plane. 
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assuming  that  M = 1.5  and  M = 2.5,  respectively.  These  results  illustrate  the 
competition  between  the  effects  of  decreasing  pumping  at  large  r,  which  tends  to 
decrease  the  intensity  as  r increases,  and  decreased  saturation  in  the  exit-ray 
region,  which  tends  to  increase  the  intensity. 

6.  Summary 

We  have  described  a simple  geometrical  model  by  which  one  can  obtain 
intensity  distributions  in  a high-Fresnel -number  unstable  resonator  containing 
saturable  gain.  This  model  is  not  a substitute  for  the  more  exact,  but  complex, 
methods  of  calculation  which  include  diffraction,  but  it  does  serve  as  a means  of 
obtaining  rough  estimates  of  intensities  in  a quick,  easy,  and  inexpensive  manner. 
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Interior  and  exit  rays  in  the  symmetric  unstable  resonator 
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Fig.  3.  Qualitative  tra 
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Fig.  7. 
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Fig.  12. a. 


2 

Confocal  resonator.  GO  = 4 exp  [-0.5(r/a  - 0.75)  ] for 
r > 0.75a.  Output-mirror-plane  intensity.  M = 1.5 
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Fig.  12. b.  Confocal  resonator.  GQD  = 4 exp  [-0.5(r/a  - 0.7 5)2]  for 
r > 0.75a.  Output-mirror-plane  intensity.  M = 2.5 
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SECTION  III 


THEORY  OF  MODES  IN  A LOADED  STRIP  CONFOCAL  UNSTABLE  RESONATOR 

Summary 

A method  is  described  for  calculating  the  modes  of  strip  confocal  unstable 
resonators,  using  a technique  that  is  practical  and  economical  even  at  high 
Fresnel  numbers.  The  present  theory  is  based  on  work  of  Horwitz.  In  contrast 
to  Horwitz' s theory,  the  present  theory  does  not  give  divergences  at  the  shadow 
boundaries.  It  also  applies  not  only  to  bare  resonators,  but  to  resonators  in 
which  there  is  arbitrary,  slowly  varying  spatial  dependence  of  the  gain.  It 
appears  straightforward  to  extend  the  present  theory  to  apply  to  three-dimen- 
sional resonators  with  rectangular-edge  mirrors.  A variety  of  numerical 
results  for  the  symmetric  modes  is  presented  graphically. 

In  recent  years  unstable  optical  resonators  have  found  wide  appli- 
cation for  high  power  gas  lasers. ^ These  resonators  offer  the  advan- 
tages of  large  mode  volume,  automatic  output  coupling,  and  good  mode 
control.  The  configuration  of  choice  in  most  applications  is  the 
positive-branch  confocal  unstable  resonator,  in  which  the  output  from 
the  resonator  is  in  the  form  of  a collimated  beam.  In  performing  theo- 
retical analyses  of  unstable  resonator  modes,  there  are  several  principal 
difficulties.  One  is  that,  because  the  light  emerges  around  the  edges 
of  at  least  one  of  the  output  mirrors,  diffraction  at  the  mirror  edges 
has  to  be  taken  into  account.  A second  factor  that  has  to  be  taken  into 
account  is  the  transverse  spatial  dependence  of  the  gain.  The  spatial 
dependence  of  the  index  of  refraction  also  is  sometimes  a significant 
factor. 
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Theoretical  determination  of  resonator  modes  gen  -ally  requires 

(21 

computer  calculations.  In  a geometrical  theory'  ' very  little  computer 
time  is  required,  but  the  results  are  crude  and  not  valid  at  low  Fresnel 
numbers.  In  theories  which  include  diffraction  the  costs  of  computing 
are  often  large  or  even  prohibitive.  Quite  elaborate  computer  codes 

based  on  the  fast  Fourier  transform  (FFT)  algorithm  have  been  devel- 

(ql 

oped.'J/  These  codes  account  for  diffraction,  gain,  and  index  of 
refraction  effects,  but  are  not  practical  to  use  at  high  Fresnel  numbers 

\ 

(Fg^  t 10).  Even  at  low  Fresnel  numbers  it  is  expensive  to  obtain 
results  for  a sizable  sample  of  resonator  parameters.  For  a bare 
resonator  several  methods  are  available  besides  the  FFT  method.  These 
include  Fox-Li  type  iterative  calculations^  and  matrix  techniques.^ 
In  the  matrix  techniques  the  non-Hermi tean  integral  equation  describing 
the  bare  resonator  is  approximated  by  a matrix  eigenvalue  equation,  and 
the  eigenvalues  and  eigenfunctions  are  calculated  numerically.  The 
matrix  techniques  give  the  higher-loss  modes  in  the  resonator,  whereas 
the  FFT  and  Fox-Li  methods  ordinarily  give  only  the  lowest-loss  mode. 
With  these  latter  methods  extra  work  is  required  to  determine  the  other 
modes,  and  it  is  not  clear  how  to  do  this  for  loaded  resonators.  None 
of  these  three  methods  works  well  at  large  F^. 

In  this  paper  we  describe  a new  method  of  calculating  modes  of  a 
strip  confocal  unstable  resonator;  i.e.,  a 2-dinensional  resonator  with 
cylindrical  mirrors.  It  is  well-known  that  the  integral  equation 
describing  a 3-dimensional  bare  resonator  with  rectangular-edged  mirrors 
can  be  factored  into  two  equations  for  strip  resonators.  While  this 

i 
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factorization  ordinarily  is  not  valid  for  loaded  resonators,  it  nevertheless 
appears  straightforward  to  generalize  the  present  theory  to  such  3-dimensional 
resonators.  Our  theory  requires  little  computer  time.  It  works  efficiently  for 

in  the  range  from  about  one  to  several  hundred,  provided  that  the  magnifica- 
tion M is  not  too  close  to  one.  It  yields  the  high-loss  modes  as  well  as  the 
lowest-loss  mode.  Our  theory  is  an  extension  of  a theory  of  Horwitz.^  In 
contrast  to  Horwitz's  theory,  our  theory  does  not  give  divergences  at  the  shadow 
boundaries.  It  also  applies  not  only  to  bare  resonators,  but  to  resonators  in 
which  there  is  arbitrary  slowly  varying  spatial  dependence  of  the  gain.  However, 
at  the  present  stage  of  development  it  does  not  include  index-of-refraction 
effects. 

1 . Basic  Formalism 

The  geometry  of  the  confocal  unstable  resonator  is  shown  in  Figure  13.  It 
consists  of  a concave  mirror  and  a convex  mirror  aligned  so  as  to  share  a common 
focal  point.  (Actually  it  is  a focal  line  in  3-space  if  the  mirrors  are  cylin- 
drical.) For  expository  purposes  we  take  the  point  of  view  that  the  concave 
mirror  is  at  the  left  end  of  the  cavity  and  the  convex  mirror  is  at  the  right 
end  of  the  cavity.  The  cavity  length  is  D,  and  the  focal  point  lies  a distance 
d to  the  right  of  the  convex  mirror.  The  convex  mirror  (also  called  the  output 
mirror)  has  half-size  a,  and  the  concave  mirror  is  assumed  to  be  large  enough 
to  reflect  virtually  all  of  the  light  traveling  to  the  left.  In  practical  terms 

a concave  mirror  half-size  (M+l)a  is  sufficiently  large,  where  M = (D+d)/d  is 

(2  7) 

the  resonator  magnification.  In  geometrical  theories  ' * ' the  light  traveling 
back  and  forth  in  the  resonator  is  assumed  to  consist  of  a left-traveling  cylin- 
drical wave  emanating  from  the  focal  point  and  a right-traveling  plane  wave. 

We  use  two  coordinate  systems  in  our  analysis.  One  is  a Cartesian  system. 
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This  has  a z coordinate  which  has  its  origin  at  the  concave  mirror  and  increases 
to  the  right  and  an  x coordinate  which  has  its  origin  on  the  resonator  axis  and 
increases  upwards.  The  other  is  a polar  coordinate  system  centered  at  the  focal 
point.  This  has  a radial  coordinate  p and  an  angular  coordinate  6,  such  that 
points  exactly  to  the  left  of  the  focal  point  have  9=0  and  such  that  6 
increases  in  the  clockwise  direction.  We  assume  that  the  paraxial  approximation 
is  valid,  so  that  in  the  arguments  of  slowly  varying  functions  we  may  take 

z = D + d - p , (1) 

X = 9p  . 

Sometimes  we  use  a mixture  of  the  two  coordinate  systems.  Thus,  we  take  the 
saturated  gain  G(p,x)  to  be  a slowly  varying  function  of  the  longitudinal  coor- 
dinate p and  transverse  coordinate  x. 

We  consider  G(p,x)  to  be  a given  function  (up  to  an  additive 
function  of  p alone).  The  other  half  of  the  problem,  that  of  calculating 
G in  terms  of  the  intensity  distribution  in  the  laser,  is  not  treated 
here.  A complete  solution  for  the  laser  steady  state  could  be  made  by 
iteratively  calculating  the  intensity  distribution  and  the  gain  distri- 
bution in  turn  until  both  converge  to  a self-consistent  solution.  We 

(2) 

have  applied  this  iterative  approach  to  our  earlier  geometrical  theory, 
where  we  obtain  convergence  after  about  five  iterations.  Since  the  main 
objective  in  our  present  work  has  been  to  test  the  validity  of  our  approach 
in  a simple  2-dimensional  geometry  before  proceding  to  a 3-dimensional 
theory,  we  have  not  considered  it  worthwhile  to  do  an  iterative  calcula- 
tion for  the  strip  resonator.  We  assume  that  the  gain  is  symmetric. 
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G(p,x)  = G(p,-x),  so  that  the  resonator  modes  are  either  symmetric  or  anti- 
symmetric. Also,  we  assume  that  the  mirrors  are  perfectly  aligned  and  perfectly 
reflecting.  We  take  the  field  E to  obey  the  scalar  Helmholtz  equation. 

The  effective  Fresnel  number  is  given  by 

Feff  = a2k/(4Tid)  , (3) 


where  k is  the  wave  vector.  This  parameter  is  called  the  equivalent  Fresnel 

( i 345) 

number  (N  ) by  some  authors.  ’ ’ ’ ' Resonator  mode  losses  are  well-known  to 
eq 

exhibit  a quasi-periodicity  in  F such  that  single-mode  operation  is  favored 
near  half-integral  values  of  F 

We  write  the  field  in  the  cavity  in  the  form 


E = f(z,x)  exp  [ik(z+Md)]  - g (p,e)  eik“  . (4) 

The  optical  phases  match  on  the  mirror  surfaces.  If  we  use  Eq.  (4)  in  the 
Helmholtz  equation,  we  obtain  the  following  paraxial  equations  for  f and  g: 


2 i k (af/az  - Gf)  + 32f/3x2  = 0 , (5) 

2 i k (3g/3p  - Gg)  + p'232g/3e2  = G . (6) 

Also,  since  E must  be  zero  on  the  mirrors,  we  obtain  from  Eq.  (4)  the  boundary 
condi tions 

f (0,MdQ ) * (Md)~*  g(Md,G)  , (7) 


g(d,e) 


dJsf(D,de)  e2ik0  , i 0 1 < a/d 
0 , | 9 | > a/d 


(8) 
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Before  discussing  the  effects  of  diffraction,  we  examine  briefly 

(2) 

how  our  geometrical  theory  applies  to  the  strip  resonator.  Our  geo- 
metrical theory  is  obtained  by  neglecting  the  last  terms  on  the  left- 
hand  sides  of  Eqs.  (5)  and  (6).  This  gives  the  rate  equations 


3f/3z  = Gf  , 

(9) 

3g/3p  = Gg 

(10) 

Note  that  we  define  the  gain  G as  a logarithmic  derivative  of  amplitudes, 
rather  than  intensities.  The  logarithmic  derivative  of  intensities  is 
obviously  2G.  It  turns  out  that,  if  is  large,  the  effect  of 
diffraction  simply  adds  a low-amplitude  high-frequency  ripple  to  the 
geometrical  solution.  Sufficiently  close  to  the  axis  it  is  consistent 
with  the  geometrical  theory  to  regard  G,  f,  and  g as  independent  of  x 
and  0.  Setting  x = 9 = 0 and  taking  the  Wronskian  of  Eqs.  (9)  and  (10), 
one  obtains 

f(0,0)  g(Md.O)  * f (0,0)  g(d,0)  . (11) 

Using  Eqs.  (7)  and  (3)  in  Eq.  (11)  gives 

e2ikD  M'?g2(Md,0)  = g2(d,0)  . (12) 

Comparing  Eq.  (12)  with  the  solution  of  Eq.  (10)  , 

g(p,0)  * g(d,0)  exp[ [ G(p',0)  dp']  , (13) 

J d 
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for  o = Md,  we  obtain 


y - 1 , (14) 

where 

U = M" * e2ikD  exp[2  j G(p',0)  dp']  . (15) 

Jd 

Thus,  gain  distributions  which  are  consistent  with  a steady  state  in 
the  cavity  must  satisfy  Eq.  (14).  If  we  include  the  effects  of  diffrac- 
tion, Eq.  (14)  is  no  longer  valid.  We  shall  show  that  Eq.  (14)  may  be 
replaced  to  a good  approximation  by  the  condition  that  u be  the  root  of 
a polynomial. 

Let  us  define 


fMd 

Q = G(p\0)  dp' 

Jd 


(16) 


Also,  let 

G(p,x)  * Q/D  + G' (p.x)  . (17) 

We  regard  G'(p,x)  as  a given  function.  The  integrated  gain  of  G'  along 
the  axis  is  zero.  Q cannot  be  specified  beforehand,  but  must  be  calcu- 
lated from  the  theory;  Q is  known  once  u has  been  determined.  It  turns 
out  that  in  order  to  calculate  the  field  on  the  output-mirror  plane,  only 
the  integral  of  G'  along  various  rays  traversing  the  resonator  need  be 
known.  Thus,  there  are  many  gain  distributions  G'  which  give  the  same 
field  on  the  output-mirror  plane.  For  instance,  if  we  add  to  G'  any 
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function  of  p alone  whose  integral  over  the  cavity  is  zero,  then  the  integrated 
gain  is  unaffected. 

In  this  paper  we  assume  that  the  laser  operates  at  steady  state  with  a 
single  frequency.  However,  the  case  where  G'  h 0 is  mathematically  equivalent 
to  the  case  of  decaying  modes  in  a bare  resonator.  We  accordingly  refer  to  a 
laser  having  G'  = 0 as  a bare  resonator,  even  though  we  understand  that  G cannot 
be  zero  if  a steady  state  exists. 

If  we  wish  to  account  for  the  effects  of  diffraction  at  the  mirror  edges, 
it  is  not  possible  to  neglect  the  second  derivatives  in  Eqs.  (5)  and  (6).  Dif- 
fracted rays  emanate  from  the  mirror  edges  and  bounce  back  and  forth  across  the 
cavity.  The  phases  associated  with  these  rays  have  rapid  transverse  oscillations. 
If  one  looks  at  the  field  associated  with  these  rays  well  inside  the  geometrical 
shadow  boundaries,  then  the  field,  aside  from  slowly  varying  amplitudes,  consists 
of  cylindrical  waves  emanating  from  the  mirror  edges  and  the  virtual  images  of 
the  mirror  edges.  These  image  points  (pn,±xn)  are  given  by  the  equations 


dM2n  , 

(18) 

aMn  , 

(19) 

where  n takes  on  all  integral  values.  Clearly  these  points  lie  on  a parabola. 
For  n<  0 the  points  lie  to  the  right  of  the  cavity,  for  n = 0 the  points  are 
the  output  mirror  edges,  and  for  n > 0 the  points  lie  to  the  left  of  the  cavity 
(see  Fig.  14)  if  we  wish  to  calculate  E out  to  the  region  of  tne  snaaow  boun- 
daries and  beyond,  then  the  two  points  for  n = 0 are  to  be  regarded  as  edges 
of  an  aperture  illuminated  by  cylindrical  waves  coming  from  the  image  points  for 
which  n < 0.  Similarly  the  two  points  for  n = 1 are  to  be  regarded  as  edges  of 
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of  an  aperture  illuminated  by  cylindrical  waves  emanating  from  the 
image  points  for  which  n > 1. 

In  accordance  with  the  above  discussion,  we  write  the  amplitudes 
f and  g appearing  in  Eq.  (4)  in  the  form 

N 

f(z.x)  * l (Pn-p)'2  [exp{i2-(x-xn)2/(pn-p)}  fn(p,x) 
n=l 

± exp{i£(x+xn)2/(pn-p)}  fn(p.-x)]  + f(p,x)  , (20) 

-N+l 

p"^(p»0)  = l (p-Pn)’ ‘ Cexp{i|{0-xn/pn)2/(-^  + ^~)} 

n=0,-l  ,-2,. . . n 

X gn(p,x)  ± exp{i|(0  + ^)2/(~  + ^-)>  9n(p.-x)] 

+ p"  ‘ g(p.x)  . (21 ) 

It  Is  consistent  with  edge  diffraction  effects  to  take  f , f,  gn,  g to 
be  slowly  varying  functions,  whereas  f and  g have  rapid  oscillations, 
particularly  at  large  F Cur  chief  objective  is  to  calculate 
f(D,x). 

The  phases  appearing  in  Eqs.  (20)  and  (21)  are  the  Fresnel  approxi- 
mation to  the  phases  of  cylindrical  waves  emanating  from  points  (pn,txn). 
These  phases  are  rapidly  oscillating  if  | n J is  small,  but  approach 
the  constant  phase 

exp(ika2/2d)  = exp(2iriFgff ) (22) 
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as  j n j becomes  large.  The  number  N is  to  be  chosen  large  enough  so  that 
this  limit  has  essentially  been  reached  when  jn[  = N.  (Horwitz  finds  that 
taking  = 250  F ^ gives  a sufficiently  large  value  of  N.)  Thus,  we  have 
rays  which  emanate  from  the  mirror  edges  and  make  N round  trips  of  the 
cavity.  The  rays  associated  with  f and  g,  on  the  other  hand,  originate 
in  the  core,  as  do  the  rays  of  our  geometrical  theory.  Equations  (20) 
and  (21)  are  supposed  to  apply  well  inside  the  shadow  boundaries.  The 
way  in  which  these  equations  are  modified  to  apply  in  the  outer  part  of 
the  resonator  is  described  later.  The  ± signs  refer  to  the  symmetric 
and  antisymmetric  modes. 

It  is  fairly  obvious  [and  may  be  verified  by  inserting  the  terms  of 
Eqs.  (20)  and  (21)  into  Eqs.  (5)  or  (6)]  that,  provided  G is  slowly 
varying,  the  amplitudes  fn,  f,  g , g all  obey  rate  equations  of  the  form 


3fn  3g 

3s  G fn  * 3 s 


n _ r „ 
- - G gr 


etc. 


(23) 


where  3/3s  indicates  the  directional  derivative  along  the  direction  of 
ray  propagation. 

Besides  the  rate  equations,  we  need  to  specify  the  boundary  conditions 
satisfied  by  the  various  amplitudes.  The  boundary  conditions  at  the  con- 
cave mirror  are  easily  found  by  setting  p = Md  and  9 = x/(Md)  in  Eqs.  (20) 
and  (21)  and  equating  coefficients  of  like  exponentials  when  these  expressions 
for  f and  g are  used  in  Eq.  (7).  This  gives 

(Pn-Md)‘*fn(Md.x)  = (Md  - Pl_n)"^gi.n(Md,x)  , 

n * 1,2,. ...N  , (24) 
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The  boundary  conditions  at  the  right  end  of  the  cavity  require  more 
care.  The  total  field  must  vanish  on  the  mirror,  whereas  the  left 
traveling  part  of  the  field  must  vanish  for  |x|  > a.  Let  us  suppose 
that  G = 0 in  a region  just  inside  the  right  end  of  the  cavity.  This 
assumption  is  made  for  mathematical  convenience,  and  our  final  results 
probably  do  not  depend  on  it.  With  this  assumption  g(p,9)  can  be  given 
for  p near  the  right  end  of  the  cavity  in  terms  of  g(d,0): 


g(p.e)  - ^[2ikm/(l-l)]* 


de'  g(d.e')  exp[£-  (e-e')2/(l  - J)] 


(26) 


Using  Eq.  (8)  in  Eq.  (26)  and  setting  n = 9'd/a,  we  get 


g(p,9)  = 2V  [21l«r/(;"  • J)3 i ad_i:  e2lk° 


P d 


f 

-1 


x | dn  f(D,an)  exp  [^-{9  - an/d)2/(^-  - j)] 


(27) 


The  correct  sign  of  the  square  root  is  that  for  which  the  real  part  is 
positive.  We  insert  the  specific  form  of  f(D,an)  in  Eq.  (20)  and  carry 
out  the  integration  over  n using  an  approximate  formula  developed  in 
the  appendix  of  Horwitz's  paper.  For  our  purposes  this  formula  can  be 
written  in  the  form 


I 


rl 


exp  [i (An2  - <:Bn  + 0]  F(n)  dn 

J 
-1 

* ( 1 ir/A ) * F(B/A)  exp  [i(C  - B2/A)] 

(28) 

+ F(-l ) [2i (A  + B)]'1  exp  [i (A  + 2B  + C)] 

+ F(1 ) [2i (A  - B)]'1  exp  [i (A  - 2B  + C)]  . 

In  Eq.  (28)  it  is  assumed  that  t and  S are  large  and  that  B/A  lies  in  the 
domain  of  integration  and  not  close  to  n 1 ±1 . It  is  also  assumed  that 
F(n)  is  slowly  varying.  Physically,  the  first  term  on  the  right  hand 
side  of  Eq.  (28)  corresponds  to  reflection  and  the  other  two  terms 
correspond  to  diffracted  light  emitted  from  the  mirror  edges.  Equation 
(28)  becomes  singular  if  A = ±B.  These  singularities  are  at  the  positions 
of  the  geometrical  shadow  boundaries.  The  singularities  in  Eq.  (28)  can 
be  readily  eliminated,  as  we  show  later.  However,  for  now  we  restrict 
our  attention  to  points  well  inside  the  shadow  boundaries. 

Once  the  integration  over  n in  Eq.  (27)  has  been  carried  out  and 
the  form  of  g(p,6)  in  Eq.  (21)  has  been  substituted  for  the  left  hand 
side  of  Eq.  (27),  the  slowly  varying  functions  multiplying  the  various 
rapidly  oscillating  exponentials  are  equated  term  by  term.  After  a 
certain  amount  of  algebra  one  obtains  the  following  results: 

A /N 

The  functions  f , f,  gn>  g obey  the  boundary  conditions  which  one 
would  expect  on  the  basis  of  geometrical  optics. 
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(29) 


e2lkD  (pn  - d)'?  fn(d,x)  = (d  - P_n)"f  g_n(d,x)  , 

n * 1 N-l  , 

o i i/  n -1  9 * 

e {(pN  - d)  * exp  (ika  /2d)  [fN(d,x)  ± fN(d,-x)]  + f(d,x)} 

(30) 

= d"^g(d,x) 

The  last  two  terms  on  the  right  hand  side  of  Eq.  (28)  lead  to  an  expression 
for  gQ(p,x),  the  amplitude  of  the  diffracted  wave  on  its  first  trip  across 
the  cavity.  This  expression  is  only  correct  near  the  output  mirror  and  must 
be  multiplied  by  the  exponential  of  the  gain  integral  according  to  Eq.  (23) 
In  order  to  be  valid  over  the  whole  length  of  the  cavity.  One  obtains 

N 

g.(p,a+(p-d)a)  = -e2^kD  (i/2mk)'  {£  (p  - d)_* 

0 n=l  n 

x [exp  (i|(xn-a)2/(pn-d)}  fn(d,a)  (an^  - a)'1 

(31) 

± exp  (i|(a+xn)2/(pn-d)}  fn(d,-a)  (an^  - a)-1] 

+ f(d,a)  (a/d  - a)  H exp  [ G(p*  ,a+(p' -d)a)  dp'] 

Jd 

In  this  equation  a is  the  angle  with  respect  to  the  axis  of  the  ray  from 

the  mirror  edg<-  at  (d,a)  to  the  point  where  a is  evaluated.  Thus,  if 

o 

the  ray  hits  the  concave  mirror  at  position  x^,  then  a = (xf.-a)/[d(M-l )]. 

The  angles  a ^ are  the  angles  with  respect  to  the  axis  of  rays  emanating 


45 


from  the  points  (p  n,±x  ) and  striking  the  output  mirror  edge  at  (d.a).  Speci- 
fically, a ^ = ( a/d ) ( 1 * M n)  \ These  are  just  the  angles  of  the  shadow  boun- 
daries. These  angles  rapidly  approach  a/d  as  n becomes  large  if  M is  not  close 
to  one.  These  shadow  boundaries  are  illustrated  in  figure  2. 

Equation  (31)  describes  the  propagation  of  the  amplitude  of  the  diffracted 
wave  on  its  first  trip  across  the  cavity.  By  solving  Eq.  (23)  and  applying  the 
boundary  conditions  (24)  and  (29)  the  diffracted  wave  may  be  propagated  through 
N round  trips  across  the  cavity.  The  solutions  of  Eq.  (23)  are 


fn(Md-z,x-(x-xn)z/(yd-pn))  = fn(Md,x) 


(32) 


x exp  [ dz'  G(Md-z‘ ,x-(x-xn)z'/(Md-pn))]  , n = 1,...,N  , 
■'o 


9n(p.x+(p-d)(x-xn)/(d-pn))  = gp(d,x) 

(33) 

x e p [[  G(p' ,x+(p’-d)(x-x  )/(d-p  ) ) dp']  , n = -1 ,-2 , . . . ,-N+l  . 
Jd  n n 


If  we  were  to  make  an  iterative  calculation  in  which  G had  to  be  determined 

in  terms  of  the  intensities  in  the  cavity,  tnen  we  would  need  to  evaluate 

f_  and  g„  throughout  the  cavity  (on  a suitable  grid).  However,  for  the 
r.  r» 

present  purposes  we  concentrate  our  attention  on  the  amplitudes  f (d,x), 
f(d,x)  from  which  the  right-traveling  field  of  the  laser  can  be  determined 
on  the  output  mirror  plane,  according  to  Eqs.  (4),  (20)  and  (21).  If  we 
set  z = 0 in  Eq.  (32)  and  p = Md  in  Eq.  (33)  and  apply  the  boundary 
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conditions  (24)  and  (29),  then  fn(d,x)  can  be  expressed  in  terms  of 
f^(d,x)  by  the  equation 


fn(d,x)  = e2i(n-1)kD  M(n-1)/2  Cr,|(x)/r1(sn(x))]  f1(d,sn(x))  , (34) 


where 


s„(x) 


1-M'2  ,.1-n  . 

M x + o a 


1-M 


-2n 


1-M 


-2n 


(35) 


In  Eq.  (34)  we  have  introduced  the  function  Pn(x),  which  is  defined  as  the 
exponential  of  the  integrated  gain  along  a ray  which  leaves  the  edge  of  the 
output  mirror  at  point  (d,a),  makes  n round  trips  across  the  cavity  and 
hits  the  output  mirror  plane  at  point  (d,x).  Explicitly  Tn(x)  is  given 
by  the  formula 


rn(x) 


n fMd 
exp  { l [ 
m=l  • 


Ml-m 

G(p'  > oT  (-M~nx  + a) 


1-M 


-2n 


+ £lM 


, u-l-n+m 


(x-M~na,)  dp'  + 


1-M 


Md 


■n+m 


M M r. 

G(p' , 327T  (X_M"  a) 


1-M 


+ ^ ^-7-  (-M'nx  + a))  dp']  } . 

0 1-M‘^n 


(36) 


The  first  integral  in  Eq.  (161  gives  the  integrated  gain  on  the  m'th  trip 
to  the  left;  the  second  integral  gives  the  integrated  gain  on  the  m'th 
trip  to  the  right.  If  we  examine  the  ray  paths  g’ven  in  Eq.  (36),  we 
see  that  those  rays  which  remain  in  the  resonator  for  many  bounces  enter 
and  leave  the  core  along  paths  nearly  antiparallel  or  parallel  to  the 
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rays  of  the  geometrical  theory.  All  except  the  first  and  last  few  trips 
take  place  close  to  the  axis. 

We  define  r(i'(x)  in  the  same  way  as  Tn(x),  except  that  G'  is  to  be 
substituted  for  G in  Eq.  (36).  Then 

rn(x)  = e2nQrn'(x)  . (37) 


Also  note  that 

V - If*  e2Q+2ikD 


(38) 


The  function  T '(x)  converge  to  a limit  function  as  n becomes  large, 
n 

since  rays  very  close  to  the  axis  do  not  contribute  significantly  to 
the  gain  integrals.  This  convergence  is  the  principal  reason  for 
defining  G'  as  we  do. 

We  can  relate  f-j(d,x)  to  fp(d,±a)  by  starting  with  Eq.  (31)  and 
making  one  round  trip  of  the  cavity.  This  gives 


ll 

f,  (d,x)  = -e2lkD  M-  r,(x)  (i/2irk)‘  { £ d"^(M2m-l)'‘ 
1 1 m=l 


r - rika2  l-M-n\  r fJ  , ,a  1 

x [exp  I-5-7-  f (d.a)  m 

2d  t m d ^ 


xM_1-aM“2  v-1 


1+M 


1+M 


d(l-M_z) 


)' 


iexp'{1#  rS}  fm(d*-a)  1 


1-M 


d 1-M'm 


xM~] -aM~2  N-1 
d(l-M-2) 


)“] 


(39) 


+ f(d,a)  (|  . )-l}  . 

d d(l-M  z) 
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Inspection  of  Eq.  (39)  reveals  that  the  major  part  of  the  x dependence  of 
f-|(d,x)  is  of  the  form  r.|(x)/(a  - x/M).  Accordingly , if  we  define  q(x)  by 


f 1 ( d » x ) = e2lkD  M*  rT(x)  q(x)/(a-x/M) 


(40) 


then  q(x)  will  be  nearly  constant.  Setting  x = ±a  in  Eq.  (34),  substituting 
into  the  righthand  side  of  Eq.  (39),  and  using  Eq.  (40),  we  obtain 


q(x) 


-(id/2*ka2)*  { l (M2m-1 )”  * [exp  1^} 

^•u  , , ,.-m 


N 

V 

L 

m=l 


1+M 


X e2'""  M"'2  LW  .) 

m 1+M'm 

, (i+M~m)  (1  - )-' 

(1+M~m) (a-x/M) 

e2'«  .f'2  r (.a)  q(^  a) 


2d  1-M"m 

x (l-M'"1)  (1 aM-rcQ+M  2)  j-l-j 

(l-M-m) (a-x/M) 


1-M 


-m 


(41) 


+ d*a?(d,a)  (1-M-2)  } . 


T ♦ 4 r 

• w • J llwUv  J JUi  J 


mCv;  to  say  iGors  about  th6 


functions  ? and 


y- 


If  we 


start  with  Eq.  (25),  propagate  ? and  § to  the  right  end  of  the  cavity 
using  the  rate  equation  (23),  and  substitute  into  Eq.  (30),  we  obtain 

A 

the  following  functional  equation  for  f(d,x): 


<(PN-d) 


-l 

A. 


1ki  /2d  [fN(d.x)  ± fN(d,-x)]  + f(d,x)} 
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(42) 


2ikD  -A  . , 
x ei1KJ  M ‘ exp  { 


Md 


G( p * ,Mx)  dp'  + 


Md 


G(p\xp'/d)  dp'} 


= f(d,Mx) 


From  Eq.  (35)  we  see  that  s^(x)  = a/M  for  large  N.  Therefore,  from 
Eqs.  (34)  and  (40)  we  obtain 

fN(d,x)  = e2i(N'1)kD  M(ri-1)/2  CrN(x)/r7 (a/M)]  f^d.a/M) 

(43) 

- e2lNkd  MN/2  rN(x)  q(a/M)/[a(l-M"2)]  . 


Substituting  Eq.  (43)  into  Eq.  (42)  and  using  Eqs.  (37)  and  (38)  gives 
for  large  N 

(d'<  exp  (ika2/2d)  uN  [rN*(x)  ± rN* (-x)]  q(a/M)  a'1  (l-M-2)'1  + f(d/x)} 


x u exp  { 


Md 


G'(p'  ,Mx)  dp'  + 


Md 


G(p' ,xp'/d)  dp'} 


= f(d,Mx)  . 


(44) 


We  choose  N sufficiently  large  that  r„'(x)  has  reached  its  limiting  form. 

il 

N 

This  may  or  may  not  require  a larger  N than  the  Horwitz  criterion  M = 
250  ?effi  depending  on  the  particular  problem  being  solved.  By  a con- 
sideration of  the  ray  paths  for  large  N,  one  can  see  that  r^'(x)  is  a 
symmetric  function.  The  solution  of  Eq.  (44)  for  the  antisymmetric 

A 

modes  is  simply  f = 0.  The  solution  of  Eq.  (44)  for  the  symmetric  modes 
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is 


f(d,x)  « [uN+1/(l-u)]  d*  ‘ exp  (ika2/2d)  2rN'(x)  q(a/M)  a'1  (1-M‘2)-1. 

(45) 

This  may  be  verified  directly  by  substituting  Eq.  (45)  into  Eq.  (44). 

In  deriving  Eq.  (45)  we  have  taken  the  limit  N -*•  °°  in  certain  of  the 
factors  (such  as  the  phase).  In  numerical  calculations  it  is  actually 
simpler  not  to  use  the  limiting  form.  Then  f(d,x)  has  the  same  structure 

as  the  other  terms  in  Eq.  (20)  (specialized  to  z = D)  except  for  division 
by  (1-u).  In  the  subsequent  analysis  we  do  sometimes  use  the  limiting 

form  and  sometimes  do  not,  as  is  convenient. 

If  we  set  x = a in  Eq.  (45)  and  insert  into  Eq.  (41),  we  get  an  equation 


In  which  the  only  unknown  amplitude  is  q(x).  However,  further  approxima- 
tion is  necessary  in  order  to  get  a tractable  polynomial  equation  for  y. 
We  wish  to  make  an  approximation  in  which  q(x)  is  constant  and  can  be 
divided  out  of  Eq.  (41).  There  is  no  unique  way  to  do  this.  One  way 
is  to  neglect  M~m  compared  to  one  everywhere  where  M m appears  in 
Eq.  (41)  except  in  the  phases.  A second  way  is  to  approximate  the  angles 
of  all  the  shadow  boundaries  by  a/d.  With  either  of  these  approximations 
it  is  consistent  to  regard  q(x)  as  constant.  Horwitz  uses  the  second 
approximation.  For  the  most  part,  so  do  we.  However,  there  seems  to 
be  no  obvious  reason  why  one  approximation  is  better  than  the  other. 
Fortunately,  in  these  cases  in  which  we  have  done  the  calculation  both 
ways,  the  results  are  in  very  good  agreement  with  each  other. 

Combining  Eqs.  (41)  and  (45),  dividing  by  q,  using  the  second 
approximation  just  described,  and  also  using  Eqs.  (37),  (38)  and  (3), 
we  obtained  the  desired  polynomial  equation  tor  y: 
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1 


, , N+1  , 

= ' 5T  (i/2Feff)?  {J,  ^ {27Ti  Feff  V Bm'  5 

m=  l 

x rm'(a)  = exp  t2xi  Feff/Sm>  6m?  !•„'(-«)]  u"  (46) 

x (1/(1 -y)  for  symmetric  modes  if  m = N+l)}  , 

where 

8m  = (1-M'm)/(1+M 'm)  . (47) 

If  we  had  used  the  first  approximation,  the  factors  6m  2 and  sm2  would  have 
been  absent  from  Eq.  (46).  In  the  case  of  a bare  resonator  [r  '(±a)  = 1] 

Eq.  (47)  becomes  the  same  as  Horwitz's  polynomial,  aside  from  differences  due  to 
an  overall  complex  conjugation  between  our  theory  ahd  Horwitz's  theory. 

The  rays  along  which  the  gain  is  integrated  to  give  r '(±a)  are  rays  which 
leave  the  mirror  edge  at  point  (d,a),  make  n round  trips  of  the  cavity,  and  end 
up  at  point  (d,±a).  Some  of  the  ray  paths  are  shown  in  Figure  3 for  a resonator 
with  M = 2.5.  Those  rays  which  return  to  point  (d,a)  do  so  by  retracing  in 
reverse  direction  the  path  which  they  travel  on  their  first  n one-way  transits. 
Because  of  the  assumed  symmetry  of  the  gain,  the  total  integrated  gain  along 
those  ray  paths  which  arrive  at  (d,-a),  as  well  as  those  which  arrive  at  (d,a), 
is  just  twice  the  integrated  gain  along  the  first  n one-way  transits.  One  can 
see  in  Figure  15  that  the  ray  paths  for  large  n approximate  the  ray  paths  of  the 
geometrical  theory. 

Once  Eq.  (46)  has  been  solved  for  u,  the  amplitudes  fn(d,x)  and  f(d,x)  are 
easily  solved  for.  From  Eq.  (40)  we  get  (putting  q = 1) 

f 1 (d ,x ) = PMr1 ' (x)/(a-x/M)  . (48) 
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The  amplitude  f(d,x)  is  given  by  Eq.  (45)  (with  q = 1)  for  the  symmetric  modes 
and  is  zero  for  the  antisymmetric  modes. 

The  total  field  amplitude  f(D,x)  of  the  right-traveling  wave  on  the  output- 
mirror  plane  is  found  by  using  Eqs.  (49)  and  (45)  in  Eq.  (20).  However,  this 
does  not  give  a valid  solution  except  well  inside  the  shadow  boundaries.  In 

particular,  Eq.  (48)  diverges  at  x = Ma.  In  the  next  section  we  show  how  to 

(3) 

eliminate  this  divergence. 


2.  Extension  Throuih  the  Shadow  Boundaries 

So  far  we  have  restricted  our  attention  to  points  well  within  the  shadow 
boundaries.  In  order  to  extend  our  results  to  large  |x|  we  need  a better 
approximation  to  the  integral  I than  that  given  in  Eq.  (28).  We  write  I in  the 
form 


1 

I * exp[i(C-B2/A)]  exp[iA(n-B/A)2]  F(n)  dn 

J-1 


X 


(A/iir) 


exp(iAr)  d s 


(50) 


where  we  have  simply  multiplied  by  an  integral  whose  value  is  one.  We  interpret 
n and  ? as  Cartesian  coordinates,  so  that  the  integration  in  Eq.  (50)  is  over 
an  infinite  strip.  Next  we  convert  to  polar  coordinates  r and  t centered  at 
5 = 0,  n = B/A,  as  illustrated  in  Figure  16.  First  suppose  that  ! B/ A ! < 1.  Then 


I * (A/Itt)'*  exp [ i (C-B2/A) ] { 


*/2 


-it/2 


dt 


(1-B/A)  see* 


F(n)rdr  exp(iAr2) 


+ 


3tt/2  (-1-B/A)  sect 

d$  [ F(n)rdr  exp(iAr2) }. 

Jit/2  J0 


(51) 
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We  wish  to  obtain  an  approximate  expression  for  I which  is  valid  if 
F ( n ) is  slowly  varying.  We  do  this  by  integrating  by  parts  according  to 


2 

r exp(iAr  ) F(n) 


jj)r  fW  IF  leMUrZ)1 

llr  texP(iAr2)  F(n)]  - exp(iAr2)  ^7  F(n)> 


and  neglect  the  term  involving  the  derivative  of  F(n).  Then  Eq.  (51) 
becomes 

I = (1tt/A)  ^ exp[i(C-B2/A)]  f (B/A)  - (i/^A)1* 

(53) 

x [ei(A-2B+C)  p(1)  R(A(l-B/A)2)  + ei(A+2B+C)  F(-l)  R ( A ( 1 +B/A ) 2 ) ] , 


where  the  function  R(p)  is  defined  as 


rir/2  2 

R(p)  = exp(ip  tan  ;)  d$. 

-tt/2 


(54) 


The  function  R(p)  is  essentially  the  Rubinowicz  line  integral  appropriate 
(9) 

to  a straight  edge.  It  is  possible  to  derive  corrections  to  Eq.  (53) 
involving  the  first  and  higher  derivatives  of  F(n)  at  n 5 B/A,  1,  and  -1 
However,  we  assume  that  Eq.  (53)  is  an  adequate  approximation. 

If  we  now  consider  the  case  B/A  - 1 , we  obtain  instead  of  Eq.  (53) 
the  expression 


I = -(i/4rrA)‘  [-ei(A_28+C)  F(l)  R(A(1-B/A)2) 
+ ei(A+2B+C)  F(.1)  R(A(1+B/A)2)3  . 


(55) 


54 


Since  R ( 0 ) = ir,  we  see  that  Eqs.  (53)  and  (55)  give  the  same  result  when 
A = B.  An  expression  for  I which  is  valid  for  all  values  of  B/A  is  given 
by 

I = (itt/A)  ‘ exp[i(C-B2/A)j  F ( B/A)  e(A2-B2)  - (1/4* A)*5 

x [e1 (A-2B+C)  p(1)  R(A(1_b/A)2)  e(A_B)  (56) 

+ el(A+2B+C)  F(-l)  R(A(1+B/A)2)  e(A+B)]  , 


where 


0(p) 


P > 0 , 


0,  p < 0 


(57) 


e(p) 


' 1.  P > 0 , 
-1 , p < 0 . 


(58) 


Equation  (53)  can  be  reduced  to  Eq.  (28)  by  making  a stationary  phase 
approximation  of  Eq.  (54)  valid  for  large  p: 


R(p)  = (iir/p)  ‘ 


(59) 


The  physical  meaning  of  Eq.  (56)  is  that  diffraction  of  a cylindrical  wave 
by  a slit  gives  rise  to  a total  field  that  is  the  sum  of  three  contribu- 
tions: the  undiffracted  field,  chopped  off  at  the  geometrical  shadow 

boundaries,  and  two  cylindrical  waves  emanating  from  the  edges  of  the 
slit.  These  latter  waves  are  not  quite  proper  cylindrical  waves,  however, 
since  they  behave  peculiarly  near  the  shadow  boundaries.  In  fact,  they 
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are  discontinuous  ac  the  shadow  boundaries  in  just  such  a way  that  the 
total  field  is  continuous.  The  right-traveling  field  in  the  resonator 
may  be  viewed  as  a set  of  cylindrical  waves  emanating  from  the  points 
(p  ,±x  ) n = and  a plane  wave.  These  waves  are  diffracted  by 

a virtual  aperture  at  (pj,±Xj).  Similarly,  the  left-traveling  field 
may  be  viewed  as  a set  of  cylindrical  waves  emanating  from  the  points 
(pn,txn),  n = -1,...,-N+1,  and  (0,0).  These  waves  are  diffracted  at 
the  output  mirror  (pQ»ix0)-  In  accordance  with  Eq.  (56)  the  amplitudes 

A 

f 2 , f ^ and  f at  p = d are  to  be  chopped  off  at  the  shadow  boundaries. 

The  shadow  boundaries  of  the  plane  wave  are  at  x = ±M a.  The  shadow 
boundaries  of  the  cylindrical  waves  are  at  x = Ma(l  + M‘n~^)/(1  + M-n+^) 
and  x 1 -Ma(l  - M~n”^)/(1  - M"n+b.  These  differ  significantly  from 
x = ±Ma  if  n is  small.  Equation  (49)  is  altered  to  read 


fn(d’x) 


un  Tn'(x)  (a-sn(x)/M)-1  e(Ma(l+M‘n_1)/(l+M‘n+1)  - x) 
x 0(x+Ma(l-M‘n'1)/(l-M‘n+1))  , n = 2.....N  . 


(60) 


Equation  (45)  is  modified  similarly  by  multiplying  by  8(Ma-x)  6(x+Ma). 
We  see  that  Eq.  (60)  is  well-behaved  and  that  fn(d,x)  is,  in  fact, 
slowly  varying.  The  singularity  where  sn(x)  = Ma  occurs  at  x = Mna, 
which  is  well  outside  the  shadow  boundaries,  provided  that  M is  not 
close  to  one. 

It  remains  now  to  determine  f^(d,x).  To  this  end  we  calculate 
the  edge-diffraction  contributions  to  Eq.  (27),  using  Eqs.  (56)  and 
(20).  We  find  that 
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(61) 


, , 1 J 2ikD  -A  . 

90(P*X)  = - 27  d“  e { L (P-<*)'0n  ' (P-P.n)  ' 

n=  l 

X [exp  (xn-a)2/(pn-d)}  fn(d,a)  R(^  (an^-a)2 

x (p-d)(pn-d)  (p-P_n)_1  Pn_1)  c(an^+^-a) 

^ 2 

± exp'  {^-  [ — > fn(d,-a)  R(y-  (an^'^-a)2  (p-d)  (en*d)  (p-p_n ) _1  Pn_1) 

2 

X e(an^-ci)]  + (p-d)  ‘ p~t  [elka  /2d  PN~  ® (fN(d»a)  ± fr|(d,-a)) 

+ f(d,a)]  R{~-  - a)2  (p-d)  p-1)  e(|-  - a)} 

[P 

x exp' { I G(o‘ , a+(o'-d)a)  do'} 

•<d 

Equation  (61)  reduces  to  Eq.  (31)  well  inside  the  shadow  boundaries,  where 
Eq.  (59)  is  valid,  except  that,  to  be  consistent  with  our  truncation 
scheme,  we  extend  the  summation  in  Eq.  (61)  only  up  to  n = N-l  and  lump 

A 

the  term  for  which  n = N together  with  the  term  containing  f(d,a).  To 
calculate  f-|(p,x)  we  regard  the  points  (p^,±x^)  as  edges  of  a virtual 
aperture  illuminated  by  cylindrical  waves  emanating  from  points  (pn,±xn), 
n * 2,...,N  and  by  a plane  wave.  Near  p = p^  the  total  amplitude  f(z,x) 
is  given  by  the  integral 

1 i rxi 

f(Md-p,x)  = ^7  [-Zirki/fp.j-p)] i j ' dx*  f(Md-p1fx') 

~X1  (62) 
x exp  [4f-  (x-x')2/(p1-p)] 
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If  we  let  n = x'/x-j  and  use  the  form  of  f(Md-p-j  ,x^n)  implied  by  Eq.  (20) 
then  we  can  use  Eq.  (56)  to  carry  out  the  integration  over  q.  The  function 
f.|(p,x)  comes  from  the  second  term  in  Eq.  (56),  which  describes  edge 
diffraction  at  q = 1 . We  obtain 


Vp’x>  = ’ 2¥  {J2  (Pl'P)Mpn"P)  [exp  {i2  Pn-Pl 


• (xrxnr 


x ^(p-i  »x-j)  R(2  (p1-p)(pn-P1)(pn-p)  1 (pj.p  + ) 


x,-x  X 
X e(-L—  + 

prp  p, 


1 ”*n  k ^ ^ 

^)±exp’^W-}fn(pr-xi) 


(63) 


i x,-x  x,+x  2 x.-x  x.+x 
■l/l  . 1 n>  ^ 1 , 1 n, 


x R(f  (p,-p)(p  -pJ(p-p)''  (-1—  + -L-i)  ) e(-i—  + — — — ) ] 
2 1 n 1 n P1  -P  Pn-P1  P-j -p  Pn"P] 


+ (p-|-p)  ‘ f(Pl  ,x1 ) R(j  (p1  -p)  (^7^-)  ) c(x1-x)} 


If  p is  not  close  to  p^ , then  Eq.  (63)  must  be  modified  by  multiplying  the 
right-hand  side  by  the  exponential  of  the  gain  integrated  along  the  rays 
emanating  from  (p^,x^).  In  particular,  at  p = Md  one  must  multiply  by 


exp  l 


Md 


G(p\  a+(p'-d)  dT^fy)  dp') 


in  order  to  match  f^(Md,x)  to  gQ(Md,x)  on  the  concave  mirror: 


(p.|-Md)~  * f.j  (Md,x)  = (Md-d)‘?  gQ(Md,x) 


(64) 
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At  o = d the  right  hand  side  of  Eq.  (63)  is  to  be  multiplied  by  T^x). 

A 

The  coefficients  fn(p^,±x^)  and  f(p^,x^)  appearing  in  Eq.  (63)  are  found 
by  substituting  Eqs.  (63)  and  (61)  into  Eq.  (64).  One  gets 


fn(pT4xl) 


e2lkD  fn-l  (d’±a)  • n = . 


(65) 


,,..2, 


f(prx1)  = e2,kD  [elka  /2d  PN"‘(fN(d,a)  ± fN(d,-a))  + f(d,a)] 

(66) 


The  values  of  f (d,±a)  and  f(d,a)  are  gotten  by  setting  x = ±a  in  Eqs. 

(49)  and  (45): 

fn(d,±a)  = un  Mn  rn'(±a)  (1  ± M'n)/[a(l-M“2)]  , (67) 

f(d,a)  = d"‘a_1  (1-M"2)'1  exp(ika2/2d)  (liljiy  (a)uN+1/(l-u).  (68) 

If  we  use  Eq.  (57)  (for  n = N)  a :d  Eq.  (68)  in  Eq.  (56)  and  use  the 
limiting  form  of  f^  for  large  N,  we  get 


f(Pl,Xl)  * e2ik0  M’*  exp(ika2/2d)  d'V^l-M-2)"1 


* (l±l)rN'(a)uN/(l-u). 
Substituting  Eq.  (67)  into  Eq.  (65)  gives 


(69) 


fn(r»1,±x1 ) = e2ikD  un_1  Mn']  r^Ua)  (l±M'n+1  )/[a(l-M'2)]  . (70) 
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Next  we  substitute  Eqs.  (69)  and  (70)  into  Eq.  (63),  set  p = d,  and 
multiply  the  right-hand  side  of  Eq.  (63)  by  F^(x).  This  gives 


• n * 1 N+l  « i 

f,  (d,x)  = - ± r/(x)  M(l-M-Z)'?a_l  l (MZn-l)_? 

1 **  1 n=2 

x un  Mn  [exp  {i2m  Feff  } rn-l(a)  n+M_n+1) 


x R(2tt  F 


M'2-M'2n 


eff  (1-M'2n)(l-M'2) 


l+M'1'"  x^x 
(M - a}  } 


1+M 


(71) 


x e(M  - j)  - exp'  {i2T  F 


1+M 


1+M~n+1  } , (_a) 

eff  -|_M-n+l  n-1 


x ( 1 -M-n+1 ) R ( 2rr  Ffiff 


M-2-M~2n 


(1-M’2n)(l-M"2) 


(M 


1-M'n+1 


x 


e-M 


l-H'"'1 


(l/(l-j.)  for  symmetric  modes  if  n = N+l). 


The  total  amplitude  f(D,x)  is  found  by  setting  z = D in  Eq.  (20)  and 
substituting  Eqs.  (50),  (45)  multiplied  by  3 functions , and  (71)  into 
Eq.  (20).  After  some  algebraic  simplification  using  Eq.  (35)  and  assuming  x > 0 
we  obtain  the  basic  result: 


f(D,x) 


Y tu"  (l-i'l"2n)  * (1-M*2)'1 

n=2 


[exp'  { i 2tt 


r n-M^x)2, 
eff  1.M-2n 


rn'(x) 


-1 


x (1-M~nx)  9(^T  - ft  1 exp  {12tt  Feff  1.M-2n 


-n-1 


(1+M~nx)2- 


Tn'(-x) 
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„ .„-n  ,-l  /1-M 

x (1+M  x)  e( — - ») J 

l-M  1 11 

.-In  2 


- [(exp  {12tt  Feff  ^ iy(x)  (1  -M_2)_1{ Un(l-M"2n)' 


I +M 


(72) 


x R(2t  F 


1-M 


-2n+2 


(1-M  iM)(l-M*‘)  1+M" 


r1+M'1-n  x,2,  /1+M~n~^  xv 

L „J.T  - mJ  ) e(  _n+]  - M) 


1+M 


+ exp  { i 2tt  F 


1+M 


•n+1 


eff  n-1 


•}  r’  (-D  (i-M'n+1) 


1-M 


-2n+2 


1-M"n_1  x^2.  x, 


x R(2ir  F 1=5- r r - gT)  e(-^ ryy  - £)]}) 

eff  (1-M-2  ) (1-M"2)  1-M  M 1-M  A 


± (same  expression  with  x replaced  by  -x ) ] > x ( 1 / ( 1 -u ) 'For 
the  symmetric  modes  if  n = N + 1) 


In  Eq.  (72)  we  have  set  a = 1,  i.e.,  x is  measured  in  units  of  a.  It  is 
straightforward  to  check  that  Eq.  (72)  is  continuous  at  the  shadow  boundaries. 
However,  this  does  not  rule  out  the  possibility  of  significant  discontinuities 

in  slope.  It  is  possible  that  such  discontinuities  occur  at  values  of  M very 
close  to  one.  In  tne  numerical  calculations  wnicn  we  have  carried  out,  wnere 

we  take  M * 2.0  or  M * 2.9,  there  is  no  noticeable  lack  of  smoothness  in  f(D,x). 

For  the  symmetric  modes  the  coefficients  of  un  in  Eq.  (72)  become 
independent  of  n for  n sufficiently  large.  Note  that  the  last  term 
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(for  n = N + 1)  is  to  be  divided  by  1-y.  The  exact  value  of  N is 
unimportant,  so  long  as  it  is  sufficiently  large.  Further  increases 
in  N simply  add  more  high-loss  unphysical  modes  by  increasing  the  degree 
of  the  polynomial  in  Eq.  (46).  For  the  lowest-loss  mode  u is  typically 
quite  close  to  one,  especially  near  half-integral  values  of  Feff,  so  that  the 
term  which  has  (1-y)  as  a divisor  is  considerably  larger  than  the  other  terms. 

For  the  antisymmetric  modes  the  coefficients  of  un  in  Eq.  (72) 
approach  zero  as  n becomes  large.  Physically  this  means  that  waves 
emanating  from  the  two  mirror  edges  bounce  back  and  forth  and  spread 


• m f ^ 1 U 
ui  i U i ' 


■»  r.  r o 1 


by  destructive 


interfcrenc 


3.  .lumen cal  Calculation  of  the  Symmetric  Mode*; 

In  this  section  we  present  a variety  of  solutions  for  symmetric 
modes  of  strip  confocal  unstable  resonators.  Since  the  lowest  loss  mode  is 
generally  symmetric,  the  antisymmetric  modes  are  of  less  interest,  and  we  do 
not  discuss  them  further  here. 

We  present  solutions  both  for  bare  resonators  and  for  loaded 
resonators.  In  the  case  of  loaded  resonators,  we  assume  that  the  gain 
6 is  independent  of  o,  extends  the  entire  length  of  the  cavity,  and 
depends  on  x according  to 

G’(x)  = A [exp(-ox2)-l] . (73) 

2 

We  take  6 = 1.5/a  and  take 


A = 1/4  in  M/[(M-l)d] 


(74) 
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or 


A = -1/8  in  M/[(M-l)d]  . (75) 

In  the  first  case  we  have  a Gaussian  gain  distribution.  The  value  of  A is 
chosen  such  that  for  the  lowest-loss  mode,  for  which  |y(  = 1,  the  gain 

G(x)  = G'(x)  + 1/4  IT)  J Mu  | / [ ( M- 1 ) d ] (76) 

falls  to  near  zero  for  large  x. 

% 

In  the  second  case  we  have  a negative  Gaussian  gain  distribution.  In  this 
case  the  gain  G(x)  for  the  lowest-loss  mode  dips  to  about  2/3  of  its  large-x 
value  as  x + 0. 

We  do  not  claim  that  either  of  these  forms  for  G'  is  realistic.  Finding 
a realistic  form  for  G would  require  an  iterative  calculation.  Such  an  iterative 
calculation  would  be  worthwhile  for  a three-dimensional  resonator  with  a 
rectangular-edged  output  mirror,  but  is  probably  not  worthwhile  for  a strip 
resonator.  The  main  point  of  the  present  work  is  to  test  the  general  soundness 
of  our  theoretical  approach  in  a simple,  two-dimensional  geometry.  Once  this 
has  been  tested,  it  appears  straightforward  to  proceed  to  three  dimensions. 

We  have  developed  two  codes,  called  ROOTS  and  MODES,  which  have  been  used 
to  calculate  the  results  presented  here.  ROOTS  uses  the  sub-routine  CPQR 
(provided  by  Kirtland  AFB)  to  calculate  the  roots  u of  the  complex  polynomial 

in  Eq.  (46).  For  the  examples  presented  the  degree  of  the  polynomial  which  we 

use  ranges  between  9 and  13.  In  the  case  of  a loaded  resonator  it  is  necessary 
to  calculate  r '(±a).  This  is  done  by  straightforward  numerical  integration  of 

Eq.  (36).  We  do  this  by  using  parabolic  interpolation  and  Simpson's  rule  inte- 

gration. We  specify  G'(x)  on  a 25  x 25  grid  extending  the  whole  length  of  the 
resonator  and  from  x = 0 to  x = 4.8a.  We  have  also  tried  using  a 49  x 49  grid 
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and  get  essentially  the  same  results.  ROOTS  does  not  operate  at  maximum  efficiency, 
since  after  each  increment  in  F CPQR  starts  from  scratch  and  recalculates  all 
of  the  roots.  Horwitz's  technique  of  using  previously  calculated  roots  as  trial 
solutions  and  omitting  calculation  of  roots  correspondi ng  to  high-loss  modes 
would  make  ROOTS  operate  much  more  efficiently.  At  present,  using  a CDC  6400 
computer,  a calculation  of  roots  at  30  values  of  F typically  requires  about 
7.9  sec  execution  time  for  a bare  resonator  and  14.6  sec  for  a loaded  resonator. 

Part  of  the  reason  for  the  longer  time  required  for  a loaded  resonator  is  that, 
for  most  of  the  cases  we  have  considered,  a higher  degree  of  polynomial  is  re- 
quired for  a loaded  resonator  in  order  to  make  the  r's  converge. 

MODES  calculates  and  plots  the  intensities  and  phases  of  the  resonator 
modes.  For  an  M = 2.9  resonator,  for  example,  this  involves  calculation  at  401 
values  of  x/a  using  a step  size  of  0.01.  This  takes  about  10.3  sec  for  a bare 
resonator.  For  a loaded  resonator  one  must  also  calculate  rn'(x).  Since  the 
r's  are  slowly  varying,  we  calculate  them  using  a step  size  of  0.2  and  interpo- 
late onto  the  finer  grid.  It  is  inefficient  to  calculate  the  r's  using  Eq.  (36). 
Instead,  we  calculate  them  successively  using  the  equation 


rn'(x)  = exp' { 


Md 


J-n 


G'(p',  — 2^  (-M‘nx  + a) 

d 1 -M~^n 


+ £’  M'1  , „-n 


/•Md 


d 1-M‘2n 


(x-M  a))  do'  + 


G(p\ 


1-M 


-2n 


x (x-M’na)  + — — 5 - (-M’nx  + a))  do'} 

a 1-M_‘n 


(77) 


x rn  1 ' ( — [(-M1_2n  + M_1 )x  + (M1'"  - M‘n_1)a]) 
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The  time  needed  to  calculate  the  r's  is  about  5.9  seconds.  Since  the  r's  do  not 
depend  on  F ^ or  on  the  particular  mode  being  calculated,  it  is  often  useful  to 
calculate  several  different  intensity  and  phase  curves  in  one  run  while  calculat- 
ing the  r's  only  once.  The  functions  rn’(x)  are  shown  graphically  in  Figure 
for  the  case  of  a Gaussian  gain  distribution  in  a resonator  with  M = 2.9.  For 
small  n the  functions  are  asymmetric,  but  they  approach  a symmetric  limit  function 
as  n becomes  large. 

The  function  R(p)  which  is  defined  in  Eq.  (54)  must  be  calculated  in  order 
to  determine  f(D,x)  according  to  Eq.  (72).  This  function  may  be  expressed  in 
terms  of  the  Fresnel  integrals  f(x)  and  g(x)  by 

R(p)  = "(1-1)  {g[(2p/*)b]  + if  [(Zp/tt)*5]  . (78) 

The  functions  f and  g are  well  behaved  and  may  be  conveniently  evaluated  by  means 
of  a rational  approximation  given  in  the  mathematical  1 i terature. 

The  results  presented  here  are  of  high  quality  in  the  sense  that  the  grids 
used  are  sufficiently  fine  to  give  results  of  good  accuracy  and  high  resolution. 

The  amount  of  computation  time  required  is  nevertheless  not  large.  If  coarser 
grids  were  used,  the  time  required  would  be  even  smaller. 

In  Figures  13  and  19  we  show  plots  of  |X[  versus  Fpff,  where  |X|  = l/lu|.  The 
results  are  shown  in  this  form  to  facilitate  comparison  with  Horwitz's  results. 
Figure  18  is  for  a resonator  with  M = 2.9  and  15.1  < F nCe  < 11.0.  Figure  19a  is 
for  a bare  resonator.  Figure  19b  is  for  a resonator  containing  Gaussian  gain 
distribution,  and  Figure  19c  is  for  a resonator  containing  a negative  Gaussian  aain 
distribution.  Note  that  the  presence  of  a Gaussian  gain  distribution  enhances 
the  separation  of  the  lowest-loss  mode  from  the  other  modes,  while  a negative 
Gaussian  gain  distribution  has  the  opposite  effect. 


65 


Figures  20  throuqh  22  show  intensity  and  phase  distributions  |f(D.X)l  and 
arg(f(D,x)]  in  the  output-mirror  plane  for  a resonator  with  M = 2.9  and  F ^ = 16.4. 
The  tick  marks  on  the  positive  side  of  the  horizontal  axis  indicate  the  positions 
of  the  shadow  boundaries.  Figure  20  shows  the  lowest-loss  mode  of  a bare  resonator. 
Figure  21  shows  the  next-to-lowest-1 oss  eymmetric  mode  of  a bare  resonator.  The 
intensities  in  these  two  figures  are  in  complete  agreement  with  Horwitz's 
Figure  24  in  the  range  0<  x<  a where  comparison  is  possible.  Note  the  large 
peak  in  intensity  in  Figure  21  near  the  shadow  boundaries  the  small  intensity 
near  the  axis.  Also  note  that  while  the  phase  in  Figure  20  is  nearly  constant  out 
to  the  shadow  boundaries,  the  phase  in  Figure  21  fluctuates  considerably.  The 
apparent  (but  unreal)  discontinuities  in  phase  are  merely  the  result  of  requiring 
the  phase  to  lie  between  -180°  and  +180°,  normalized  with  respect  to  the  phase 
at  the  center  of  the  mirror.  Intensities  are  normalized  so  that  the  largest 
intensity  in  the  range  0 < x < a is  unity. 

Figure  22  is  for  the  lowest-loss  mode  of  a resonator  containing  a Gaussian 
gain  distribution.  This  may  well  be  the  first  calculation  yet  performed  for  a 
resonator  containing  distributed  gain  and  having  this  high  a Fresnel  number.  As 
one  might  expect  with  this  sort  of  gain  distribution,  the  intensity  dies  off  as 
x increases  well  before  x gets  near  the  shadow  boundaries.  However,  the  fine 
structures  of  the  intensity  distributions  in  Figures  20  and  22  are  quite  similar. 

The  phase  distributions  in  Figures  20  and  22  are  also  nearly  the  same. 

Figures  23  through  25  show  intensity  and  phase  distributions  for  a resonator 
with  M = 2.0  and  Feff  = 8-4'  These  figures  are  respectively  for  the  lowest-loss 
modes  of  a bare  resonator,  a resonator  with  a Gaussian  gain  distribution,  and 
a resonator  with  a negative  Gaussian  gain  distribution.  These  results  are  for  a 
value  of  F ^ where  FFT  calculations  are  feasible,  and  it  would  be  of  considerable 
interest  to  have  results  of  FFT  calculations  to  compare  with  our  results. 
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Unfortunately,  we  have  not  found  anything  comparable  in  the  published  literature. 

We  have  carried  out  a number  of  other  mode  calculations  which  space  does 
not  permit  us  to  present  here.  These  calculations  show  no  obvious  breakdown  of 
the  theory  at  small  ?eff  Our  calculations  have  used  values  as  low  as  F ^ = 2.9. 
Horwitz  claims  good  agreement  with  other  methods  of  calculation  down  to  F ^ = 1. 
None  of  these  results  of  other  workers  extend  to  the  domain  x > a or  to  loaded 
resonators.  We  have  not  tried  values  of  M smaller  than  2.0.  Our  theory  breaks 
down  for  M too  close  to  one,  at  least  as  a practical  matter,  since  the  degree  of 
the  polynomial  required  approaches  infinity  as  M -*■  1 . However,  it  should  be 
possible  to  get  quite  a bit  closer  to  M = 1 before  this  limitation  on  the  theory 
becomes  important. 

An  extension  of  the  present  theory  to  three-dimensional  resonators  with 
rectangular  mirrors  appears  straightforward.  Guided  by  the  fact  that  the  paraxial 
wave  equation  for  a bare  resonator  can  be  factored,  we  propose  replacing  Eqs.  (20) 
and  (21)  by  equations  containing  double  sums,  one  sum  for  each  transverse  dimen- 
sion. The  waves  associated  with  the  terms  in  these  double  sums  would  not  simply 
be  cylindrical  waves,  but  would  be  waves  having  one  caustic  on  an  edge  of  the 
output  mirror  or  one  of  its  images  and  another  caustic  on  an  orthogonal  edge  of 
a different  image.  The  diagonal  terms  would  correspond  to  spherical  waves  ema- 
nating from  the  corners  of  the  mirror  images. 

It  is  not  clear  whether  an  analogous  theory  can  be  developed  for  resonators 
with  round  mirrors.  The  main  difficulty  with  round  mirrors  is  that  the  edge- 
diffraction  wave  adds  constructi vely  to  give  a caustic  on  the  resonator  axis. 

The  behavior  at  large  F ^ of  resonators  with  round  mirrors  is  not  well  understood, 
since  it  has  not  been  feasible  to  carry  out  mode  calculations  for  these  resona- 
tors. No  evidence  of  the  separation  of  the  lowest-loss  mode  from  the  other  modes 
has  been  found,  but  it  is  not  known  whether  interleaving  of  the  low-loss  modes 


We  are  even  further  from  understanding  how  the 


persists  to  very  large  Fg^. 
spatial  distribution  of  gain  affects  the  mode  behavior. 
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fig.  22.  Phase  and  intensity  of  the  lowest-loss  mode  of  a resonator  containing  a Gaussian  qain 
distribution,  for  Feff  = 16.4,  M = 2.9,  y = 0.97999  - 0.0028037i 
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Fig.  24.  Phase  and  intensity  of  the  lowest-loss  mode  of  a resonator  containing  a Gaussian 
gain  distribution,  for  Fg^  = 8.4,  M = 2.0,  g = 0.9716  - 0 . 003881 5i 
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Fig.  25.  Phase  and  intensity  of  the  lowest-loss  mode  of  a resonator  containing  a negative  Gaussian 
gain  distribution,  for  F - 8.4,  M = 2.0,  p = 0.93869  - 0. 00721241 
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